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ABSTRACT  OF  THE  DISSERTATION 


Observations  of  seafloor  ambient  noise  with 
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In  1987  an  array  of  ocean  bottom  seismometers  (OBS)  was  deployed  in  the  deep 
ocean  off  southern  California  for  the  purpose  of  recording  ocean  bottom  ambient  noise.  By 
lowering  the  instruments  with  a  wire,  it  was  possible  for  the  first  time  to  create  an  array  which 
allowed  coherent  processing  of  short  wavelength  energy  in  the  frequency  range  0.05  to  30  Hz. 

The  array  consisted  of  9  working  instruments  in  place  for  30  days,  allowing  for  the 
recording  of  ambient  noise  during  several  changes  in  the  surface  environment.  Sea  surface  con¬ 
ditions  were  determined  from  the  U.S.  Navy’s  Global  Spectral  Ocean  Wave  Model  (GSOWM). 
Noise  levels  at  the  ocean  floor  have  been  correlated  with  the  swell  height  and  wind  speed 
predicted  by  the  GSOWM.  However,  the  theoretical  transfer  function  between  surface  wave 
height  and  pressure  fluctuations  on  the  bottom  did  not  predict  the  recorded  pressure  signal  to  a 
high  degree  of  accuracy. 

The  array  achieved  a  uniform  sampling  of  distances  between  sensors,  and  a  good 
estimate  of  spatial  coherence  has  been  made  up  to  156  m.  A  characteristic  pattern  of  coherence 
length  dependence  on  frequency  is  observed  throughout  the  deployment  Spatial  coherence  is 
related  to  the  wavelength  and  beam  width  of  the  energy  from  a  partially  coherent  source. 
Through  modeling  of  seismic  energy  it  has  been  possible  to  estimate  the  range  to  the  source  of 
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the  noise  energy.  The  coherence  field  is  consistent  with  a  wave  field  propagating  horizontally 
through  the  sedimentary  layer  from  a  source  at  the  edge  of  the  sedimentary  basin,  suggesting 
that  the  conversion  of  acoustic  energy  in  the  water  to  seismic  energy  in  the  ocean  floor  occurs  at 
the  rough  boundary  at  the  edge  of  the  sedimentary  basin. 

The  wavenumber  decomposition  or  beamforming  of  the  ambient  noise  presented 
some  difficulties.  Higher  mode  and  thus  faster  seismic  waves  were  predominant,  and  therefore 
the  array  did  not  have  a  sufficiently  large  aperture.  Nevertheless,  the  results  of  beamforming 
are  consistent  with  the  results  of  the  modeling  of  the  spatial  coherence. 
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CHAPTER  1 

HISTORICAL  INTRODUCTION 


Measurement  of  seismic  energy  is  limited  at  low  levels  by  ambient  noise.  The 
seismic  noise  level  varies  with  frequency  and  is  strongest  in  the  measurement  band  between  1 
and  10  seconds  period.  Energy  in  this  band  is  known  as  microseisms  and  has  a  long  history  of 
study.  It  was  noted  near  the  turn  of  the  century  that  a  strong  background  signal  existed  below  1 
Hz  [Wiechert,  1904],  Over  forty  years  ago  it  was  realized  that  the  signal  was  related  to  storm 
activity  [ Deacon ,  1947],  There  are  a  large  number  of  references  on  the  subject  {Iyer  [1964]  pro¬ 
vides  a  good,  though  by  now  outdated,  reference  list),  mostly  consisting  of  qualitative  observa¬ 
tions. 

Miche  [1944]  evaluated  the  hydrodynamical  equations  governing  standing  waves  to 
second  order  and  identified  a  pressure  term  which  did  not  have  an  explicit  depth  dependence. 
Longuet-Higgins  [1950]  recognized  this  pressure  fluctuation  as  a  possible  source  of  the  microse¬ 
isms.  Further  development  by  Hasselman  [1963]  led  to  the  theory  of  microseisms  in  its 
presently  accepted  form.  When  sets  of  ocean  swell  interact  with  each  other,  there  is  quadratic 
pressure  term  which  can  propagate  to  the  deep  ocean  floor.  The  effect  is  strongest  for  com¬ 
ponents  traveling  in  nearly  opposite  directions.  Their  interaction  generates  an  acoustic  signal 
which  propagates  to  the  ocean  floor  and  excites  seismic  waves  if  the  wavenumber  is  equal  to 
that  of  a  resonant  mode  in  the  solid  ocean  floor.  The  nonlinear  interaction  generates  an  excita¬ 
tion  whose  frequency  is  double  that  of  the  ocean  swell. 

Haubrich  et  al.  [1963]  gave  quantitative  observational  support  for  the  2:1  ratio 
between  the  frequency  of  the  microseism  peak  and  the  ocean  swell  at  a  nearby  site.  Using  a 
three-component  seismometer  they  were  able  to  estimate  the  direction  of  the  source.  They  con¬ 
cluded  that  reflection  of  storm  swell  off  the  local  coast  line  was  resulting  in  oppositely  directed 
propagating  waves  which  in  turn  generated  seismic  energy. 
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In  the  1960’s,  treaties  limiting  the  testing  of  nuclear  weapons  relied  in  part  on  seismic 
measurements  as  a  means  of  monitoring  compliance.  This  spurred  the  development  of  the 
World  Wide  Standardized  Seismic  Network  (WWSSN)  and  an  effort  to  increase  the  sensitivity 
of  seismic  detection.  Project  VELA  Uniform  was  sponsored  by  the  Defense  Advanced 
Research  Projects  Agency  (DARPA)  and  combined  seismic  and  space  based  methods.  One  of 
the  important  methodologies  was  to  use  sensor  arrays  to  enhance  coherent  over  incoherent  sig¬ 
nals.  The  most  ambitious  array  was  the  Large  Aperture  Seismic  Array  (LASA)  in  Montana 
which  covered  an  area  200  km  on  a  side  with  525  seismic  stations. 

Another  useful  capability  of  seismometer  arrays  is  that  the  propagation  velocity  and 
source  direction  of  seismic  energy  can  be  measured  directly.  Haubrich  and  McCamy  [1969] 
used  LASA  to  determine  the  direction  and  distance  of  the  seismic  energy  source  at  the  microse¬ 
ism  frequency.  With  the  array  they  were  able  to  distinguish  surface  wave  modes  and  body 
waves.  They  demonstrated  that  the  fundamental  Rayleigh  mode  and  the  body  waves  are  gen¬ 
erated  at  the  coast  and  in  the  vicinity  of  storms  at  sea. 

During  this  same  time  period,  the  first  deployments  of  seismographs  on  the  ocean 
floor  were  taking  place.  Columbia  University  deployed  an  ocean  bottom  seismometer  (OBS) 
for  short  periods  in  the  Atlantic  in  1964  and  then  continuously  between  1966  and  1972  at  the 
base  of  the  continental  slope  of  Cape  Mendocino,  California.  It  had  been  expected  that  the 
ocean  environment  would  be  quieter  than  land  sites.  Schneider  et  al.  [1964]  discovered  that 
there  was  no  improvement  in  signal  to  noise  ratio  for  ocean  bottom  seismometers,  signal  and 
noise  are  both  enhanced  on  the  ocean  floor.  Noise  levels  at  frequencies  between  0.1  and  1  Hz 
are  greater  at  ocean  bottom  sites  than  on  land.  Latham  and  Sutton  [1966]  concluded  that  results 
were  consistent  with  the  hypothesis  that  low  frequency  noise  propagates  as  a  Rayleigh  wave  in 
the  oceanic  crust.  There  were  two  consequences  from  these  results.  First,  there  was  little  effort 
in  the  following  decade  at  making  teleseismic  measurements  from  oceanic  stations.  Second, 
OBS  technological  advancement  emphasized  measurement  at  frequencies  above  the  microseism 
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band. 

There  have  been  more  measurements  of  ambient  noise  in  shallow  water  environments 
than  in  deep  water,  both  for  reasons  of  expense  and  of  direct  application  to  the  military.  There 
is  strong  evidence  that  in  shallow  seas  various  forms  of  ambient  noise  propagate  as  interface 
waves  along  the  bottom  [Rauch,  1980;  Kibblewhite  and  Ewans,  1986]  In  sufficiently  shallow 
water,  the  pressure  fluctuation  due  to  the  passage  of  surface  wave  can  excite  the  bottom  directly. 

Experiments  based  on  OBS’s  in  deep  water  have  now  become  fairly  routine.  Recent 
improvements  in  sensors  have  made  possible  the  measurement  of  increasingly  lower  frequen¬ 
cies.  Webb  et  al.  [1986]  have  made  pressure  measurements  down  to  0.005  Hz.  Their  data  show 
clearly  what  some  call  the  "noise  notch"  between  0.02  and  0.1  Hz.  This  frequency  band  of  low 
noise  level  has  created  a  new  interest  in  making  teleseismic  measurements  on  the  ocean  floor 
and  has  spurred  further  research  into  the  study  of  deep  ocean  ambient  noise. 

Measurement  of  microseisms  on  land  have  suggested  that  they  propagate  largely  as 
low  order  surface  waves.  Low  order  modes  have  the  property  that  their  amplitudes  decay  away 
from  the  surface  to  which  they  are  bound.  In  contrast,  body  waves  of  interest  for  global 
seismology  have  no  such  dependence.  Sax  and  Hartenberger  [1964,  1965]  have  demonstrated 
that  increases  in  signal  to  noise  ratio  (SNR)  can  be  made  by  placing  the  sensor  at  depth  in  a 
borehole. 

To  test  the  analogy  for  the  oceanic  environment,  two  different  ocean  borehole  seismic 
systems  were  commissioned:  the  Marine  Seismic  System  (MSS)  [Adair  et  al.,  1986,  1987]  and 
the  Ocean  Sub-bottom  Seismometer  (OSS)  [Duennebier  et  al.,  1987],  These  have  each  been 
deployed  successfully  on  two  occasions  contemporaneously  with  OBS.  There  have  been  three 
important  observations  regarding  mieroseisms  from  these  deployments.  First,  overall  signal 
levels  are  lower  in  the  borehole  than  on  the  ocean  floor.  This  is  because  the  ocean  floor  is  typi¬ 
cally  covered  with  sediments  of  low  rigidity.  Conservation  of  energy  requires  that  seismic 
amplitudes  increase  when  propagating  through  a  low  rigidity  medium  [Duennebier  et  al.,  1987], 


4 


in  contrast  to  the  borehole  seismometers  which  were  clamped  to  crystalline  basement.  Second, 
after  making  the  aforementioned  amplitude  correction,  the  amplitude  of  the  microseism  peak 
near  0.16  Hz  is  similar  for  ocean  bottom  and  ocean  sub-bottom  locations.  This  suggested  to 
Adair  et  al.  [1987]  that  at  these  frequencies,  oceanic  noise  is  predominantly  crustal  Rayleigh 
waves  with  little  contribution  from  Stoneley  waves  which  are  trapped  near  the  water-sediment 
interface.  Third,  at  higher  frequencies,  the  noise  level  is  considerably  higher  at  the  ocean  floor, 
especially  on  the  horizontal  components.  Duennebier  et  al.  [1987]  interpreted  the  higher  levels 
as  trapped  energy  in  the  sedimentary  wave  guide. 

Very  recently,  there  have  been  significant  improvements  in  seismic  recording  sys¬ 
tems.  The  WWSSN  instruments  were  analog  and  therefore  had  limited  dynamic  range.  As  a 
consequence  they  were  configured  as  separate  sensors  for  the  long  period  and  short  period 
bands  on  either  side  of  the  microseism  peak.  Modem  digital  systems  are  capable  of  recording 
from  0.001  to  10  Hz  with  greatly  increased  dynamic  range  with  a  single  sensor.  There  is  an 
effort  to  enhance  the  WWSSN  with  other  networks  using  these  new  instruments.  Hardware 
improvements  notwithstanding,  the  global  seismic  network  is  still  largely  limited  by  that  fact 
that  seismic  stations  are  on  land.  The  majority  of  the  Earth’s  surface  which  is  covered  by  ocean 
remains  unsampled  and  therefore  global  models  are  biased.  There  is  an  initiative  to  find  ocean 
sites  for  these  broad  band  instruments  [Purdy  and  Dziewonski,  1988].  However  the  state  of 
understanding  of  ocean  floor  ambient  noise  is  still  limited.  Specifically,  there  is  little  agreement 
whether  the  noise  propagates  in  the  ocean  crust  as  Rayleigh  waves  over  large  distances  or 
Stoneley/Scholte  waves  over  short  distances.  This  question  has  an  important  bearing  on  how 
sensors  should  be  placed  at  oceanic  sites.  If  short  wavelength  Stoneley  waves  predominate, 
emplacement  of  the  sensors  in  a  borehole  below  the  skin  depth  of  the  interface  waves  would 
result  in  a  significant  signal  to  noise  ratio  improvement.  Oceanic  Rayleigh  waves  on  the  other 
hand  have  wavelengths  on  the  order  of  kilometers,  and  it  is  not  practical  to  put  a  sensor  far 
enough  away  from  the  boundary  to  have  a  useful  effect. 
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One  of  the  reasons  that  there  is  still  much  room  for  progress  in  the  subject  of  ocean 
ambient  seismic  noise  is  that  to  date,  most  measurements  have  been  made  with  individual  sen¬ 
sors  or  multiple  sensors  with  inexact  locations.  On  land,  important  progress  was  made  when 
seismic  array  data  were  examined.  Arrays  of  sensors  allowing  coherent  processing  of  signals 
have  until  recently  been  impossible  because  of  operational  difficulties.  Such  an  array  deploy¬ 
ment  was  undertaken  in  1987  and  the  following  chapters  describe  the  operation  of  the  array  and 
some  results  from  the  data  analysis. 

Chapter  2  describes  the  ocean  bottom  seismometer  array  and  the  new  technology  that 
made  its  deployment  possible.  Wire  deployments  of  the  instrument  packages  made  possible  the 
deployment  of  the  sensors  with  unprecedented  precision.  This  in  turn  permitted  an  extremely 
tight  array  which  enabled  coherent  processing  of  the  recorded  signals  in  the  frequency  band  of 
0.1  to  10  Hz. 

In  chapter  3,  the  tools  of  analysis  —  the  frequency  and  wavenumber  spectra  —  are 
described.  The  methods  by  which  these  are  computed  may  introduce  bias  into  the  results.  Two 
alternative  methods  were  used  for  both  the  frequency  and  the  wavenumber  spectra. 

The  subsequent  chapters  report  on  the  analysis  of  the  data.  Chapter  4  focuses  on 
measurements  made  by  the  individual  instruments  while  chapters  5  and  6  discusses  the  results 
from  the  array  as  a  whole.  In  chapter  4,  I  relate  changes  in  ambient  noise  level  with  time  to 
changes  in  the  surface  conditions  such  as  local  wind  speed,  local  swell  height  and  distant  swell 
height.  I  draw  some  conclusions  about  how  each  of  these  factors  affect  the  noise  level  at  the 
array  site. 

In  chapter  5, 1  present  the  analysis  of  the  spatial  coherence  measurements.  The  varia¬ 
tion  of  coherence  with  frequency  and  sensor  separation  places  constraints  on  the  mode  of  propa¬ 
gation  of  the  seismic  energy.  The  sedimentary  wave  guide  will  be  seen  to  have  a  dominant 
effect  on  the  seismic  noise  field. 
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Chapter  6  describes  the  results  from  wavenumber  analysis  of  the  array  data.  Direct 
estimates  of  the  seismic  velocity  of  the  ambient  noise,  through  the  wavenumber  spectrum,  have 
a  large  scatter  but  give  a  picture  of  the  mode  structure  which  is  consistent  with  that  observed 
through  the  coherence  model. 


CHAPTER  2 

EXPERIMENT  DETAILS 


INTRODUCTION 

Deployment  of  the  OBS  array  took  place  in  April  and  May  of  1987.  The  array  was 
situated  near  Deep  Sea  Drilling  Project  (DSDP)  Hole  469  at  a  depth  of  3.8  km  (Figure  2.1). 
The  area  is  a  400  m  thick  sedimentary  basin  with  minimal  topography  in  the  immediate  vicinity 
[Yeats,  Haq  et  at.,  1981]. 

Groups  from  Scripps  Institution  of  Oceanography  (SIO)  and  the  Naval  Ocean 
Research  and  Development  Activity  (NORDA)  participated  in  the  experiment  which  was  code- 
named  "CIRCUS".  A  total  of  thirteen  OBS  were  available:  six  from  SIO,  six  from  NORDA, 
and  an  additional  one  from  Oregon  State  University  through  NORDA.  The  instruments  were 
deployed  from  the  R/V  Melville  in  the  period  of  April  10  through  April  22.  They  remained  on 
the  bottom  for  a  month  and  were  recovered  on  two  separate  trips  with  the  M/V  Jamie  G  and  the 
R/V  Sproull. 

At  two  times  during  the  array  deployment,  mid-water  "Swallow”  floats  were  deployed 
from  the  USNS  Narragansett  for  periods  of  24  hours  in  a  related  experiment.  One  of  the  aims 
of  this  experiment  was  to  characterize  the  difference  in  noise  level  between  the  ocean  bottom 
and  at  several  depths  in  the  water  column  [D’ Spain  and Hodgkiss ,  1988]. 

OBS  ARRAY  DESIGN 

The  OBS  array  was  designed  to  enable  wavenumber  decomposition  of  the  ambient 
noise  field  in  the  spatial  frequency  range  of  extremely  slowly  traveling  interface  waves.  The 
array  therefore  had  a  small  aperture  of  150  m  and  a  wavenumber  Nyquist  frequency  of  50  km-1. 
The  wavenumber  resolution  of  an  array  depends  on  the  aperture,  or  greatest  separation  between 
sensors.  The  Nyquist  frequency  in  contrast,  depends  on  the  minimum  spacing  between  sensors. 
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These  two  conflicting  criteria  must  be  satisfied  with  a  limited  number  of  sensors. 


A  convenient  way  of  analyzing  array  performance  is  through  the  co-array.  This  is  the 
autocorrelation  of  array  element  positions  and  is  the  set  of  points  generated  by  all  the  vector 
offsets  between  each  pair  of  sensor  positions  in  the  array  [Haubrich,  1968].  The  frequency 
response  of  an  array  is  given  by 


R( k)=  £  tf2*i<k*-ki/) 
ij- 1 

When  rewritten  in  terms  of  the  co-array,  this  takes  the  form 


3.1 


M 


R(k)  =  1  +  £cos(27tk-xm). 


3.2 


m  =  1 


The  sum  is  over  all  the  nonzero  co-array  points. 

An  important  class  of  arrays  are  those  with  a  uniform  co-array.  Since  the  separation 
between  sensors  is  constant  in  this  case,  the  response  may  be  written  (for  a  linear  array) 


M 

R(k)  =  1  +  £  cos(2ji km  Ax ) 

m  =  l 


3.3 


which  takes  the  form  of  a  Fourier  series.  A  given  function  can  be  approximated  better  by  a 
greater  number  of  terms  in  a  Fourier  series,  so  it  is  desirable  to  maximize  the  number  of  points 
in  the  co-array.  It  is  not  necessary  to  have  a  uniform  array  to  generate  a  uniform  co-array,  and 
in  fact,  equivalent  co-arrays  can  be  generated  with  fewer  sensors  with  unequal  spacing.  Figure 
2.2  shows  a  simple  example  for  a  linear  array  [LaCoss  et  al.,  1968], 


The  maximum  possible  number  of  co-array  points  given  N  sensors,  is  M=N(N-\). 
This  maximum  is  achieved  for  a  minimally  redundant  array,  one  in  which  none  of  the  vector 
offsets  are  duplicated.  The  cost  of  using  a  minimally  redundant  array  is  that  the  stability  of  the 
wavenumber  spectral  estimate  is  reduced  [ LaCoss  et  al.,  1969],  It  is  possible,  however,  to  make 
up  for  the  lost  stability  by  increasing  the  stability  of  temporal  spectral  estimate  used  in  generat¬ 
ing  the  wavenumber  estimate.  Figure  2.3  shows  the  minimum  redundancy  design  for  the  array 
used  in  the  CIRCUS  experiment  along  with  its  co-array. 
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ARRAY  DEPLOYMENT 

A  compact  array  was  desired  because  of  the  high  wavenumber  of  the  signals  under 
observation.  Small  arrays  are  also  advantageous  because  the  inhomogeneity  of  the  medium  is 
reduced.  The  deployment  of  such  a  compact  array  represented  a  significant  technical  effort. 
The  small  dimension  of  the  array  required  a  placement  accuracy  of  1  to  2  meters.  Such  an  accu¬ 
racy  in  position  can  not  be  achieved  by  the  conventional  method  of  deploying  the  instruments 
by  free-fall  through  nearly  4  km  of  water.  The  solution  was  to  use  technology  from  the  "Deep 
Tow"  system  of  the  Marine  Physical  Laboratory  (MPL)  at  Scripps,  a  system  which  combines 
wire-lowercd-vehicles  and  transponder  guidance  [Spiess  and  Lonsdale  1982], 

The  deployment  of  the  CIRCUS  array  made  use  of  two  newly  developed  tow  vehi¬ 
cles,  the  "Thruster"  and  the  "Rock-scan"  sonar.  The  Thruster  has  a  propulsion  system  which 
permits  partial  maneuverability  independent  of  the  ship.  The  CIRCUS  expedition  was  one  of 
the  first  field  tests  of  this  vehicle.  Because  of  several  breakdowns,  we  were  not  able  to  use  it  for 
every  OBS  deployment.  The  backup  vehicle,  Rock-scan  sonar  provided  no  propulsion  system 
but  included  upward-  and  downward-looking  sonar  and  a  transponder. 

A  release  hook  was  developed  which  could  be  opened  on  command  from  the  ship  lab. 
The  OBS  capsules  were  attached  via  this  hook  to  a  25  m  length  of  rope  suspended  from  the  tow 
fish.  Figure  2.4  depicts  the  elements  of  the  system  (not  to  scale).  The  rope  served  as  a  shock 
absorber  during  the  lowering  of  the  OBS  and  the  tow  vehicle  into  the  water.  The  OBS  were 
most  conveniently  put  in  the  water  with  the  side  crane  while  the  tow  fish  was  deployed  and 
lowered  through  the  stem  A-frame.  Once  the  tow  fish  was  in  the  water  with  the  OBS  hanging 
beneath  it  the  system  could  be  lowered  at  a  speed  of  30  meters  per  minute.  In  3800  m  of  water 
the  trip  to  the  bottom  took  about  3  hours. 

The  feature  that  makes  the  Deep  Tow  system  an  effective  tool  for  the  emplacement  of 
instruments  on  the  bottom  is  the  creation  of  a  fixed  reference  system  with  an  array  of  acoustic 
transponders.  Acoustic  pings  are  exchanged  between  either  the  ship  or  the  tow  fish  and  the 
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transponders  which  are  moored  to  the  ocean  bottom.  The  transponders  are  interrogated  at  three 
frequencies:  10,  10.5  and  11  kHz  and  reply  at  single  frequency  of  12  kHz.  The  round  trip 
acoustic  travel  times  are  displayed  on  a  chart  recorder  in  the  ship  lab.  Three  fixed  transponders 
were  used  for  the  present  experiment.  The  desired  sensor  positions  were  marked  on  the  chart  as 
a  distance  from  each  of  the  transponders.  The  hook  was  released  when  the  proper  combination 
of  distances  was  achieved. 

Figure  2.5  shows  the  geometry  of  the  whole  system  with  all  distances  drawn  to  scale. 
Note  that  the  size  of  the  array  is  on  the  order  of  the  length  of  the  ship,  and  that  the  ship  and  the 
array  are  separated  from  each  other  by  50  times  that  length.  For  the  tow  fish  without  propulsion 
capability  the  position  of  the  bottom  of  the  wire  can  only  be  controlled  by  positioning  the  ship. 
The  R/V  Melville  is  propelled  by  cycloids  and  can  be  maneuvered  very  precisely  in  any  direc¬ 
tion.  However,  the  connection  to  the  package  at  the  end  of  the  wire  is  flexible  so  maneuvers  at 
the  surface  take  20  minutes  to  propagate  down  to  the  bottom.  Given  enough  time  and  several 
passes,  good  positioning  could  be  achieved.  Typically  it  required  2  to  3  hours  of  maneuvering 
in  addition  to  the  6  hours  round  trip  from  the  surface.  With  the  Thruster  unit,  the  bottom  time 
was  reduced  to  under  an  hour. 

The  tow  fish,  transponder  and  ship  positions  were  recalculated  using  an  iterative  least 
squares  method  [ Lowenstein ,  1969].  One  source  of  error  which  has  not  been  estimated  is  the 
amount  of  lateral  drift  by  the  OBS  capsule  after  its  release  from  the  hook.  They  were  typically 
allowed  to  fall  freely  for  the  last  50  m  to  ensure  that  the  sensor  package  would  not  be  dragged 
along  the  bottom. 

SENSOR  PACKAGE 

The  12  SIO  and  NORDA  seismographs  are  of  the  same  design  with  only  minor  differ¬ 
ences.  The  OBS  design  has  been  in  use  at  Scripps  since  1981  [Moore  et  al.,  1981].  Their 
operation  is  controlled  by  a  microprocessor  for  maximum  flexibility  of  the  recording  parame¬ 


ters. 
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The  sensors  are  a  triaxial  arrangement  of  seismometers  and  a  pressure  sensor.  The 
seismometers  are  Mark  Products  model  L-4C-3D  with  a  resonant  frequency  of  1  Hz.  Two  dif¬ 
ferent  types  of  pressure  transducers  were  used,  an  Ocean  and  Atmospheric  Systems  model  E- 
2DP  crystal  hydrophone  and  a  differential  pressure  sensor  of  a  design  by  Cox  et  al.  [1984],  The 
differential  hydrophone  has  much  better  performance  at  frequencies  below  1  Hz.  They  were 
used  on  5  of  the  OBS  in  this  experiment. 

The  sensor  output  is  fed  through  a  variable  gain  preamp  with  a  54  dB  range.  The 
amplification  level  is  determined  by  a  long  term  average  of  the  ambient  noise  independently  for 
each  channel.  This  site  turned  out  to  be  remarkably  quiet,  and  the  gain  was  usually  at  its  max¬ 
imum  setting.  The  signal  is  further  conditioned  by  a  pre-whitening  filter  and  anti-alias  filters 
with  a  cut-off  frequency  of  30  Hz.  Figure  2.6  shows  the  transfer  function  for  the  seismometer 
channels  and  figure  2.7  for  the  two  kinds  of  hydrophones. 

A  12-bit  A/D  operating  at  128  samples  per  second  converts  the  signal  to  digital  form. 
The  data  are  buffered  and  written  to  the  tape  recorder  after  the  end  of  the  recording  window  so 
that  the  mechanical  movement  of  the  tape  recorder  does  not  affect  the  signal.  The  buffer 
memory  has  a  capacity  of  128  kbyte,  or  227  seconds  for  four  channels  at  128  Hz.  There  is  a 
capacity  for  event  triggering  that  was  disabled  for  this  experiment  to  avoid  using  up  recording 
capacity.  The  data  are  written  to  tape  on  four  serial  tracks  with  8  additional  words  after  every 
64  for  time  keeping  and  error  detection.  The  normal  recorders  have  a  capacity  of  5  Mbyte  per 
channel  of  data  or  roughly  6  hours  of  continuous  data.  Five  of  the  capsules  were  deployed  with 
a  data  recorder  with  a  storage  capacity  of  20  Mbyte  per  channel. 

DATA  QUALITY 

The  deployed  positions  of  the  sensors  are  compared  to  the  desired  positions  in  Figure 
2.3.  The  realized  positions  were  within  5  m  of  the  planned  positions  with  two  exceptions.  The 
OSU  instrument  was  of  a  different  configuration  and  did  not  fit  onto  the  hook.  It  was  deployed 
by  free-fall  and  was  sited  just  to  the  west  of  the  array  center.  The  outlying  instrument  was  also 
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deployed  by  free-fall.  An  impending  shooting  schedule  did  not  leave  enough  time  for  a  wire 
deployment.  As  it  turned  out,  having  a  sensor  at  the  larger  distance  proved  quite  useful. 

Of  the  13  instruments  deployed,  9  returned  usable  data.  These  are  marked  by  the 
filled  triangles.  Of  the  remaining,  one  capsule  failed  to  resurface,  and  the  three  others  suffered 
tape  recorder  problems.  One  of  the  four  failed  instruments  was  the  OSU  instrument  and  three 
were  SIO/NORDA  instruments  which  were  set  up  with  the  experimental  high  capacity  data 
recorders.  Furthermore,  of  the  remaining  3  Cox-Webb  sensors,  the  preamps  on  2  of  them  were 
incorrectly  biased  and  yielded  poor  quality  pressure  data. 

The  actual  co-array  for  the  functioning  sensors  (Figure  2.8)  is  somewhat  irregular 
compared  to  the  desired  co-array,  but  mildly  so.  The  wavenumber  response  (Figure  2.9)  has 
higher  side  lobe  energy  than  the  ideal  one  and  a  skewed  central  peak.  The  position  of  the  outly¬ 
ing  OBS  1  was  not  included  in  the  calculation  of  these  displays.  Incorporating  that  one  would 
have  the  effect  extending  the  co-array  to  the  NE  and  SW  and  would  narrow  the  central  peak  of 
the  response  in  that  same  direction. 

It  is  reasonable  to  ask  whether  we  might  have  have  done  as  well  with  a  random  array 
deployed  inexpensively  from  the  sea  surface.  The  extreme  small  size  of  the  array  means  that  a 
fluctuation  in  the  ship  position  of  one  ship  length  would  result  in  an  sensor  position  as  far  away 
as  OBS  1.  A  random  array  of  this  sort  would  almost  certainly  have  a  much  lower  aliasing  fre¬ 
quency  and  would  not  be  capable  of  resolving  Stoneley  wave  velocities..  The  large  effort 
required  to  deploy  instruments  precisely  was  worthwhile. 

Figure  2.10  shows  the  data  coverage  for  the  good  instruments.  The  suffix  (L)  after  the 
OBS  name  indicates  the  use  of  a  low  frequency  Cox-Webb  hydrophone.  OBS  5  and  6  were  also 
fitted  with  high  capacity  recorders  and  failed  early.  We  suspect  that  the  recorders  placed  a  high 
power  burden  on  the  system. 

The  4-track  OBS  FM  tapes  are  read  by  a  computer  in  the  lab  and  stored  on  9-track 
digital  tape.  Some  of  the  records  could  not  be  processed  because  of  synchronization  errors  on 
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the  tape.  Of  the  functioning  instruments  81%  of  the  records  are  actually  available  for  analysis. 
All  in  all,  about  51%  of  the  possible  data  was  actually  recovered.  Though  one  could  hope  for  a 
better  proportion,  this  leaves  a  substantial  data  set  of  144  records  recorded  by  up  to  9  instru¬ 
ments  on  4  components. 

Aside  from  the  quality  of  the  individual  time  series,  array  analysis  requires  precise 
relative  timing  control.  The  clock  on  each  of  the  OBS  operates  independently  and  is  subject  to 
drift.  The  one  month  deployment  was  a  relatively  long  one  for  these  instruments  so  the  timing 
problem  was  serious.  The  clock  error  at  recovery  can  be  measured  against  a  reference  time,  and 
a  linear  interpolation  of  the  correction  can  be  made  for  all  recording  times.  This  was  done  and 
found  to  be  inadequate.  The  data  set  was  scanned  for  transients  appearing  on  all  or  most  of  the 
instruments  and  a  time  correction  was  made  to  align  them.  Linear  interpolation  of  the  time 
correction  between  these  infrequent  occurrences  of  transients  was  more  satisfactory.  To  use  a 
source  in  the  medium  for  this  puipose  requires  an  assumption  about  the  propagation  velocity  of 
the  signal.  These  waveforms  did  not  appear  to  be  dispersed  so  we  assumed  that  these  signals 
traveled  at  the  velocity  of  sound  in  water  and  thus  traversed  the  150  m  extent  of  the  array  in  0.1 
s.  Lacking  knowledge  of  the  source  direction  we  assumed  that  the  propagation  across  the  array 
was  instantaneous.  This  is  acceptable  if  the  array  is  to  resolve  velocities  well  below  1500  m  s_l. 
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Figure  2.1.  Map  of  the  CIRCUS  array  site.  The  array  position  is  32°  37’  N,  120°  34’  W  at  a 
depth  of  3800  m.  The  isobath  units  are  meters. 
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Figure  2.2.  Comparison  of  co-arrays  for  redundant  and  nonredundant  linear  array  designs. 
Even  though  fewer  sensors  arc  used  in  the  nonredundant  design,  the  same  lag  space  coverage  is 
achieved.  However  the  redundant  design  provides  repeated  lags. 
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Figure  2.3.  Minimally  redundant  design  of  the  CIRCUS  array  with  12  sensors. 
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Figure  2.4.  Elements  of  the  deployment  system  (not  to  scale).  The  ocean  bottom  seismometer 
capsule  is  attached  to  an  electrically  operated  release  hook  and  suspended  beneath  the  "Thrus¬ 
ter"  tow  fish.  The  Thruster  is  lowered  from  the  ship  with  0.680"  wire. 
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Figure  2.5.  Deployment  geometry  to  scale.  The  surface  ship  is  70  m  in  length,  the  array  has  an 
aperture  of  80  m.  The  array  is  at  a  distance  of  3800  m  below  the  ship. 
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Seismometer  Channel  Response 


Frequency  (Hz) 


Figure  2.6.  Response  of  the  seismometer  channel  to  ground  acceleration.  An  anti-aliasing  filter 
cuts  of  the  energy  above  30  Hz. 
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Hydrophone  Channel  Response 


Frequency  ( Hz  ) 


Figure  2.7.  Response  of  the  hydrophone  channel  to  pressure  for  both  the  commercial  crystal 
hydrophone  and  the  Cox- Webb  differential  pressure  sensor. 
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Figure  2.8.  Actual  array  configuration.  The  crosses  indicate  the  12  element  design  and  the  tri¬ 
angles  indicate  the  estimated  positions  of  the  OBS  capsules.  In  most  cases,  the  realized  position 
was  within  3  m  of  the  target.  In  two  cases  the  realized  position  does  not  correspond  to  any  tar¬ 
get  position.  The  filled  triangles  indicate  the  OBS’s  that  returned  data. 
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Figure  2.9.  Comparison  of  the  co-array  pattern  for  the  designed  array  and  the  actual  array.  The 
extent  of  the  actual  array  is  about  the  same  as  planned,  but  the  distribution  of  co-array  lags  is 


more  uneven. 


Power 
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Planned  Array  Response 


Realized  Array  Response 


Figure  2.10.  Wavenumber  response  for  the  planned  array  with  12  elements  and  the  actual  array 
with  9  working  sensors.  The  planned  array  response  has  sidelobes  of  0.2  times  the  peak  power. 
The  realized  array  has  sidelobes  of  up  to  0.5  times  the  peak  power. 
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Figure  2.11.  Data  coverage  in  time.  All  instruments  were  scheduled  to  begin  recording  at  day 
112  00:00  (CUT).  OBS  1  was  put  in  place  24  hours  late.  OBS  5  and  OBS  6  recorded  at  15 
minute  intervals  until  day  1 14  06:00  then  reverted  to  a  3  hour  recording  interval.  All  the  other 
instruments  recorded  at  6  hour  intervals.  The  capsule  names  marked  with  (L)  were  configured 
with  the  long  period  pressure  sensors. 


CHAPTER  3 

DATA  ANALYSIS  METHODS 


The  CIRCUS  array  deployment  collected  a  large  volume  of  time  series  data.  Ambient 
noise  is  essentially  a  random  process;  Haubrich  [1965]  showed  that  microseisms  are  stationary 
over  a  period  of  hours  and  approximately  normally  distributed.  The  present  data  have  not  been 
tested  explicitly  for  these  properties,  but  they  arc  assumed  to  have  similar  statistical  behavior. 
Second  order  statistical  parameters  such  as  spectra  are  appropriate  for  data  of  this  type.  The 
fundamental  tools  for  examining  the  data  collected  in  this  study  are  the  auto-spectrum  and  the 
cross-spectrum. 

POWER  SPECTRA 

The  spectral  power  density,  being  a  statistical  quantity  of  a  random  variable,  requires 
an  expectation 

P(f)  =  E[\y(f)\2)  (3.1) 

where  y  (/)  is  the  Fourier  transform  of  the  data.  The  expectation  operation  is  usually  performed 

by  some  kind  of  averaging,  of  which  a  variety  of  methods  have  been  developed.  The  different 

methods  vary  in  their  effectiveness  and  the  type  of  data  for  which  they  are  suitable.  In  this 

study,  spectra  have  been  computed  with  two  methods;  the  segment-averaged  Fourier  transform 

[Welch,  1967]  and  the  multiple  window  Fourier  transform  [ Thomson ,  1982],  The  algorithms 

are  well  known,  particularly  for  the  first  method.  They  will  be  described  here  to  emphasize  the 

difference  in  results  that  may  be  achieved  using  either  one. 

In  the  first  method,  the  time  scries  is  divided  into  many  (20  or  more)  segments  which 
may  or  may  not  partially  overlap.  Each  segment  is  tapered  with  a  suitable  function  (e.g.  a 
Kaiser-Bessel  window)  to  reduce  bias,  and  then  a  fast  Fourier  transform  (FFT)  is  performed  on 
the  segment.  The  transformed  segments  are  then  averaged  together  to  yield  the  spectral  csti- 
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mate.  If  the  time  series  x  is  sampled  at  intervals  of  A7\  then  let  xn  =x(.nAT).  The  segment 
averaged  estimate  of  the  power  spectrum  is 

P(f)  =  —Xy-Vly-V).  (3-2) 

N  «=i 

where  the  time  series  is  divided  into  N  segments  and  y*  is  the  weighted  Fourier  transform  of  the 
n-th  segment 

?,(/)  =  ?,(*  IM  A7-)  =  jf^'£wnxm+L(m-l)e-2*iim,M ,  (3.3) 

«=0 

where  M  is  the  length  of  the  segment,  L  is  the  increment  between  each  segment,  w  is  the  win¬ 
dow  function,  and  the  superscript  denotes  complex  conjugation.  The  averaging  of  segments 
is  done  to  make  the  estimate  consistent  and  to  reduce  the  variance.  There  are  two  drawbacks  to 
the  method:  one  is  that  the  resolution  possible  by  performing  the  Fourier  transform  on  segments 
of  length  M  is  much  poorer  than  would  be  possible  by  transforming  the  original  series  of  length 
(N-\)L+M .  The  other  is  that  averaging  increases  the  bias  of  the  estimate.  The  segment  averag¬ 
ing  method  is  most  effective  when  the  time  series  is  long  enough  so  that  even  when  split  into 
segments,  the  resolution  is  acceptable.  The  bandwidth  of  the  estimate  is  proportional  to  the 
reciprocal  of  the  section  length  T.  The  total  length  of  the  series  is  [(W-1)L/A/+1]7\  The  quan¬ 
tity  (1/7H(N-1)L/A/+1]7  =  ( N-\)L/M+\  may  be  called  the  time  bandwidth  product,  and  reduces 
to  N  when  L  =M .  The  time  bandwidth  product  is  a  measure  of  the  trade-off  between  the  desired 
frequency  resolution  and  the  amount  of  averaging  required  for  a  stable  estimate.  Confidence 
intervals  may  be  assigned  to  the  spectral  estimate  by  assuming  that  the  data  are  normally  distri¬ 
buted.  The  error  distribution  is  then  given  by  a  Xv  variable  with  v=2 N  degrees  of  freedom 
where  N  is  the  number  of  segments  used.  Figure  3.1  shows  some  examples  of  the  trade-off  of 
variance  and  resolution  for  some  spectra  computed  from  the  CIRCUS  data. 

The  multiple  window  method  (also  referred  to  as  the  multiple  taper  method)  is  to  be 
preferred  when  the  time  series  is  not  long  enough  to  be  divided  into  enough  segments.  The 
method  is  considerably  more  complicated  and  expensive  of  computing  resources.  It  is  based  on 
the  observation  by  Slepian  and  Pollack  [1961]  and  Slepian ,  [1978]  that  for  a  given  bandwidth 


(-WMV),  of  all  possible  time  limited  functions  of  unit  energy,  the  prolate  spheroidal  wavefunc- 
tions  minimize  the  energy  outside  the  frequency  bandwidth.  The  discrete  prolate  spheroidal 
wavefunctions  (DPSW)  are  a  natural  consequence  of  mapping  a  continuous  stationary  process 


onto  a  finite  sample  using  the  Dirichlet  kernel,  KD(f0J)= 


sinA/jc  (f0-f) 
sinrc(/0-/) 


The  DPSW, 


Uk(N,W\f  o)  arc  solutions  of  the  eigenvalue  problem 


j  Kd  (f0J  )Uk  (N  ,W-J  o)  =  K  (N  ,W)Uk  (N  ,WJ0). 


(3.4) 


The  eigenvalues  Xk  are  real,  distinct,  and  limited  between  0  and  1.  The  first  2NW  eigenvalues 
are  nearly  equal  to  1  and  then  they  decay  towards  0.  The  Fourier  transforms  of  the  DPSW  are 
called  the  discrete  prolate  spheroidal  sequences  (DPSS)  and  are  computed  from  a  discrete 
eigenvalue  problem  analogous  to  Equation  (3.4).  Slepian  [1978]  showed  that  the  estimate  using 
the  zero-order  DPSS,  v0 


(3.5) 


is  the  spectral  estimate  with  the  minimum  spectral  leakage  for  any  choice  of  window. 

Thomson's  [1982]  contribution  was  to  average  together  estimates  using  several  low 
order  sequences.  The  first  step  is  to  select  a  frequency  bandwidth  (-W,W)  or  equivalently  a 
time  bandwidth  AW,  where  N  is  the  length  of  the  data  sequence.  The  time  bandwidth  AW  con¬ 
trols  variance  and  the  resolution  of  the  spectral  estimate  and  also  the  number  of  prolate 
spheroidal  sequences  that  will  be  used  to  compute  the  estimate.  The  data  are  multiplied  with 
each  DPSS  in  turn  and  then  Fourier  transformed  yielding  K  eigenspectra,  yk .  The  spectral  esti¬ 
mate  is  given  by  an  average  of  the  eigenspectra.  An  unweighted  average  is  not  suitable  because 
the  higher  order  sequences  have  poorer  bias  properties.  Instead,  a  weighted  average  is  used 


EWly*(/)l2 


(3.6) 


where  the  weight  coefficients  are 
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and  are  computed  adaptively  from  the  data  by  an  iterative  solution  of 


(3.7) 


K-ih(Sk(f)-h^Sk(f)) 

y - ; - - - —  =  0,  (3.8) 

A  KSk(fhBk(f))2 

where  Sk  is  the  current  solution  of  the  iteration,  Bk  is  an  estimate  of  the  broad  band  bias,  and  A 
is  a  constant  used  to  normalize  for  the  case  of  white  noise.  This  algorithm  is  a  modification  of 
Thomson's  [1982]  by  Chave  [1988];  the  latter’s  notation  is  used  here.  The  equivalent  degrees 
of  freedom  may  be  given  by 

v(f)  =  2-£lkdk2(f)  (3.9) 

*=o 

which  yields  a  maximum  of  2 K  in  the  absence  of  broad  band  bias.  Thomson  [1982]  also  noted 
that  this  formulation  is  equally  suitable  for  line  spectra  as  well  as  for  continous  spectra.  A 
regression  estimate  for  the  amplitude  of  a  line  component  is 

I  Uk(N,W,0)yk(f) 

ti(/)=-^ - •  (3-10) 

I  Vk2(N,W-  0) 

An  F-test  may  be  performed  on  this  value  to  determine  if  a  line  is  significant.  This  feature  will 
be  used  briefly  in  this  study  in  the  observation  of  ship  signals. 

The  uncertainty  in  the  estimate  can  be  determined  by  a  jack-knife  procedure  [Efron, 
1982],  The  spectral  estimate  is  recomputed  with  each  successive  sequence  omitted  in  turn.  The 
difference  between  the  various  estimates  is  accumulated  and  then  converted  into  a  95% 
significance  bound.  The  resolution  and  the  variance  of  the  estimates  are  dependent  on  the 
bandwidth  chosen.  Figure  3.2  shows  estimates  of  the  noise  spectrum  made  by  the  multiple  win¬ 
dow  method  for  time  bandwidths  of  2,  4,  and  6.  Subsequent  analyses  have  used  a  time 
bandwidth  of  4,  which  constitutes  a  compromise  between  stability,  resolution,  and  computa¬ 


tional  effort. 
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By  comparing  Figures  3.1  and  3.2,  it  can  be  seen  if  high  resolution  is  not  critical,  the 
segment  averaging  method  can  yield  better  variance.  In  this  study,  when  low  frequencies  are  of 
interest,  the  multiple  window  method  will  be  used,  but  when  stability  of  the  estimate  is  desired, 
as  in  beamforming,  the  segment  averaging  method  with  a  high  time  bandwidth  product  will  be 
used. 

The  cross-spectrum  between  two  time  series 

Sij(f)  =  E[yli)(f)y^)m  (/■)],  (3.11) 

can  be  computed  directly  from  their  individual  spectra,  whether  by  the  segment-averaged 

method  or  the  multiple  window  method.  The  segment  averaged  method  requires  accumulating 

the  conjugate  product  of  the  Fourier  transforms  for  the  two  series 

=  (3.12) 

"  »=1 

The  cross-spectrum  must  be  computed  individually  for  each  segment  and  then  averaged,  rather 
than  in  the  reverse  order.  Otherwise,  the  number  of  independent  quantities  is  reduced  to  one 
(c.f.  Claerbout,  1976).  For  the  multiple  window  method,  the  cross-spectra  are  computed  from 
the  weights  and  eigcnspcctra  determined  from  the  individual  series 

Sij  (/")«£  ek(i)(f  )ep'(f  )y^{f  )  (3.13) 

k  =0 

where, 

ek(n)=dk(n)l^dJ(n)2. 

7-1 

Complex  cross-spectra  are  usefully  represented  in  polar  form:  magnitude  and  phase.  The  nor¬ 
malized  magnitude  is  the  coherence,  which  is  the  correlation  coefficient  at  each  frequency 
between  two  samples  from  and  any  jointly  stationary  processes  [ Priestley ,  1981].  It  indicates 
the  degree  to  which  the  processes  arc  linearly  independent.  The  cross-phase  is  a  measure  of  the 
propagation  delay  of  a  signal  measured  at  two  sites. 

Uncertainty  estimates  can  be  made  for  the  coherence.  For  the  segment-averaged 
method,  the  error  involves  a  complicated  combination  of  hypergcometric  functions  [ Priestley , 
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1981].  Often  the  statistic  is  reduced  to  a  threshold  determining  whether  the  coherence  estimate 
is  significantly  different  from  zero,  in  which  case  the  threshold  T=3/N ,  where  N  is  the  number 
of  degrees  of  freedom  in  the  estimate.  For  the  multiple  window  method,  it  is  again  possible  to 
form  a  jack-knife  estimate  of  the  error  by  making  estimates  of  the  cross-spectra  from  different 
combinations  of  the  eigenspectra. 

For  the  general  multi-channel  case,  the  spectral  matrix  is  given  by 

S(f)  =  E[y(f)yH  (/)]  (3.14) 

where  the  elements  of  the  complex  vector  y  are  the  spectral  estimates  at  frequency  /  for  each 

channel,  and  the  superscript  H  means  conjugate  transpose  of  the  complex  vector.  For  the  seg¬ 
ment  averaging  method,  it  takes  the  form 

S</)=^-SyI0f)y"0f).  (3.15) 

™  1*1 

where  yi  is  a  vector  of  the  transform  of  the  i-th  windowed  data  segment.  For  the  multiple  win¬ 
dow  method, 


S0f)  =  ^IE*0f)yt0f)yf  (f)Ek(f),  (3.16) 

K  *=i 

where  E*  is  a  K  xK  matrix  with  the  weights  of  the  k-th  estimate  on  the  diagonal,  and  y*  is  a  vec¬ 


tor  of  the  k-th  eigenspectra.  The  coherences  can  be  derived  from  the  spectral  matrix  by 


*</(/)  = 


Su(f)Sjj(f)’ 


(3.17) 


where  5(/  is  the  i  J-th  element  of  the  spectral  matrix.  The  spectral  matrix  is  an  estimate  of  the 
spatial  covariance  sampled  at  each  of  the  sensor  separations  and  as  such  is  also  a  fundamental 
quantity  for  the  computation  of  wavenumber  spectra,  as  will  be  seen  next. 


BEAMFORMING 

As  the  power  spectrum  in  time  yields  a  useful  statistic,  so  does  the  power  spectrum  in 
the  space  dimensions.  The  operation  of  forming  the  wavenumber  spectrum,  often  called  beam- 
forming,  is  a  spatial  decomposition  of  the  measured  wavefield  used  to  determine  the  bearing  of 
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noise-generating  sources  and  the  velocity  or  slowness  of  the  energy  propagating  from  them. 
This  can  be  seen  by  writing  the  two-dimensional  spatial  Fourier  transform 

*(x,/)  =  -M  xy(k,co)e,(k*_<“)dkdco.  (3.18) 

8tt 

The  component  corresponds  to  a  unit  amplitude  plane  wave  propagating  with  a  velocity 

c  =  lkl/co  and  bearing  (MarT1/^ /kx ,  so  the  time  and  space  domain  wavefield  is  a  weighted  sum¬ 
mation  of  all  such  components. 

From  this  point  on,  assume  that  the  temporal  transform  of  the  data  has  been  per¬ 
formed.  Unlike  Fourier  transforms  of  time  sampled  data,  spatial  data  are  typically  very  sparsely 
sampled  and  FFT  methods  are  generally  not  used.  Instead  a  sum  over  the  small  number  of  com¬ 
ponents  is  computed.  Writing  this  explicitly, 

y(\)=y(\)^8(\-\n),  (3.19) 

n  =1 

where  y  is  the  actual  and  wavefield  and  y  is  the  sampled  wavefield.  Taking  the  Fourier 
transform  of  both  sides  and  by  the  convolution  theorem, 

m  =  YQ k)  *  -ijL  £s ix-xje^dx  (3.20) 

n=sl 

y(k)  =  y‘(k)  *  (3.21) 

4Jt  B=1 

If  we  define  the  complex  array  response 

fi(k)=-rTiyk’’.  (3-22) 

4rr  b=1 

then 

Y(k)  =  Y(k)*B(k).  (3.23) 

The  measured  wavenumber  spectrum  is  the  convolution  of  the  true  wavenumber  spectrum  with 

the  array  response.  More  commonly  the  power  is  calculated  since  the  phase  of  a  two- 

dimensional  field  is  difficult  to  interpret, 

/’(k)  =  I  K(k)  I2  =  IF(k)*fl(k)l2.  (3.24) 

There  are  a  variety  of  specialized  algorithms  in  common  usage:  the  conventional 
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(Bartlett)  method  [Lacoss  et  al.,  1969],  the  maximum  likelihood  (ML)  method  [Capon,  1969], 
the  linear  prediction  method  [ Lang  and  McClellan,  1976;  Makhoul,  1982],  and  eigenvector 
methods  [Schmidt,  1986].  Analysis  of  the  CIRCUS  noise  data  have  been  performed  only  with 
the  conventional  and  ML  methods  because  they  are  relatively  independent  of  model  assump¬ 
tions. 

The  conventional  estimate  is  simply  the  modulus  squared  of  the  complex  spatial 
transform.  It  can  be  rewritten  concisely  using  matrix  notation.  Let  e(k)  be  the  "steering"  vector 
whose  elements  are  ej=e‘kx',  and  let  y(f)  be  the  vector  of  spectral  estimates  at  each  sensor. 
Then  the  discrete  spatial  Fourier  transform  of  y  becomes 


Y(f  ,k)  =  e(k)w  y(f). 

(3.25) 

and 

Pen,  if  ,k)  =  E  [  1  Y (f  ,k)  1 2]  =  ew  £  [yy*  ]e. 

(3.26) 

Substitute  the  definition  for  the  spectral  matrix  (Equation  3.14)  and 

P  conv  (f  ,k)  =  ew  Se. 

(3.27) 

From  Equation  (3.15)  it  is  evident  that  the  resolution  of  the  conventional  wavenumber  estimate 
is  limited  by  the  width  of  the  central  lobe  of  the  array  response  fi2(k).  Sometimes  a  weighted 
steering  vector  is  used  where  ej=wje,kx',  to  reduce  the  sidelobes  of  the  array  response,  but 
weighting  by  conventional  windows  cannot  decrease  the  width  of  the  central  lobe. 

Non-linear  estimates  of  the  wavenumber  spectrum  can  overcome  this  limitation. 
Capon  [1969]  derived  a  method  which  achieves  better  resolution  by  making  the  weights,  or 
steering  vectors  adapt  to  the  data.  The  weights  are  computed  for  an  optimal  window  for  each 
wavenumber,  frequency  and  spatial  signal  and  noise  structure  to  reject  both  independent  sensor 
noise  and  interfering  propagating  signals.  A  filter  is  sought  that  will  pass  a  target  plane  wave 
undistorted  with  unity  amplitude  and  also  minimize  the  total  power  output  of  the  array.  Let 
a  =  col[au  ■  ■  ■  ,aN ]  be  the  vector  of  weights.  The  first  condition  results  in  the  constraint  equation 


awe  =  1, 


(3.28) 
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and  the  second  condition  is 


p  Capon  (k)  =  min  E 

i(k) 


awSaj  . 


Use  the  method  of  Lagrangian  multipliers  to  minimize 
F  =  aw  Sa  -  pCa^  e  -  1). 

Minimize  F  by  taking  the  gradient  with  respect  to  a  and  setting  it  to  zero. 

V,  aw  Sa  -  p(aw  e  -  l)j  =0, 

2Sa-pe  =  0, 
or, 

a  =  -— S-1e. 

2 

Substitute  a  into  the  constraint  equation  (Equation  3.26)  to  solve  for  the  multiplier, 

-fK-- 

and 


(3.29) 


(3.30) 


(3.31) 


and  so 


a  = 


S~‘e 

e^S^'e 


(3.32) 


Substituting  the  weights  a  into  the  power  output  of  the  array  (Equation  3.24), 


P Capon  (f  JO  = 


r  > 

H 

f  > 

S 1e 

S 

S le 

ewS-‘e 

v  J 

e"  S”’e 

k  J 

Pcapon(fM)=-jr^r-  (3-33) 

e  S  e 

This  estimator  was  named  by  Capon  [1969],  the  maximum  likelihood  estimator.  It  does  not 
actually  correspond  to  the  estimator  derived  from  the  maximum  likelihood  analysis  but  it  is 
very  similar  in  form  [Johnson,  1982;  Kay  and Marple,  1982],  The  method  will  be  referred  to  as 
Capon’s  method  henceforth. 


The  performance  of  any  array  processing  method  can  be  gauged  by  its  resolution,  bias 
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and  stability.  These  quantities  are  sensitive  to  errors  such  as  uncertain  sensor  positions,  correla¬ 
tion  between  signal  and  noise,  and  inaccurate  assumptions  of  the  signal  model.  These  errors  are 
difficult  to  estimate  for  the  present  study.  Great  expense  and  effort  was  spent  making  the  sensor 
positions  accurate  to  within  3%  of  the  aperture.  Since  the  ambient  noise  is  the  object  of  study, 
the  distinction  between  signal  and  noise  becomes  fuzzy.  Clearly  the  signal  to  noise  ratio  is  one, 
and  there  is  a  high  degree  of  correlation  between  signal  and  noise.  The  signal  model  needs  to 
be  as  unrestrictive  as  possible. 

The  conventional  method  has  resolution  limited  by  the  beam  pattern  inherent  in  the 
array  geometry.  Lacoss  et  al.,  [1969]  showed  that  the  stability  of  the  conventional  estimate 
depends  on  the  redundancy  in  the  co-array  and  the  time-bandwidth  product  of  the  cross¬ 
spectrum  estimate,  with  an  array  geometry  factor.  The  only  model  assumption  made  is  that  the 
signal  of  interest  is  composed  of  plane  waves. 

Capon’s  method  can  achieve  improved  resolution  because  it  actually  makes  use  of  the 
coherence  structure  of  the  noise.  Capon  et  al.,  [1973]  showed  that  the  method  was  highly  effec¬ 
tive  at  removing  an  interfering  signal  that  was  highly  correlated.  But  the  question  of  correlation 
between  signal  and  noise  has  not  been  answered. 

Other  methods  can  achieve  even  better  resolution  than  Capon’s  method,  but  they  usu¬ 
ally  make  stringent  assumptions  about  the  signal  which  may  not  be  accurate  and  cannot  be 
checked.  For  instance,  the  linear  prediction  methods  assume  an  all-pole  model  for  the  signal 
and  that  the  number  of  poles  is  known.  The  eigenvector  methods  are  most  effective  when  there 
are  a  small  number  of  discrete  sources,  where  small  means  much  fewer  than  the  number  of  sen¬ 
sors. 

Figure  3.4  compares  a  wavenumber  spectrum  computed  with  the  conventional 
method  and  Capon’s  method  for  a  synthetic  signal.  The  time  series  is  a  2  Hz  sine  wave  with  a 
phase  delay  corresponding  to  a  propagation  velocity  of  0.5  km  s-1  (2  s  km-1 )  towards  the  south 
and  10%  additive  white  noise.  The  conventional  method  obviously  has  significant  sidelobe 


35 


energy  and  a  broad  main  peak.  However  the  power  level  is  approximately  correct.  The  estimate 
with  Capon’s  method  yields  a  much  sharper  spike,  but  the  peak  power  level  bears  no  simple 
relationship  to  the  input.  Wavenumber  spectra  have  been  computed  with  Capon’s  method  in 
this  study  as  a  general  rule.  However  it  is  has  proven  a  useful  check  to  compare  the  results  with 
the  conventional  method. 
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Figure  3.1.  Microseism  spectral  estimates  using  the  segment  averaged  Fourier  transform 
method  showing  the  tradc-olf  between  variance  and  resolution.  The  dashed  lines  indicate  the 
95%  confidence  interval.  The  time-bandwidth  product  reduces  to  the  number  of  segments  used 
in  the  averaging.  Because  the  record  is  of  finite  length,  more  averaging  means  less  resolution. 
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Figure  3.2.  Microseism  spectral  estimate  using  the  multiple  taper  method.  The  dashed  lines 
indicate  the  95%  confidence  interval.  The  length  of  the  Fourier  transform  is  the  same  for  each 
case  but  the  time-bandwidth  product  relates  to  the  number  of  tapers  used  in  the  calculation. 
More  tapers  results  in  a  broader  peak. 
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Separation  =  1 56  m 


IQ-1  10°  101 


Frequency  (Hz) 


Figure  3.3.  Cross-coherence  between  data  from  two  vertical  component  seismometers 
separated  by  a  distance  of  156  m.  The  dashed  lines  indicate  the  95%  confidence  interval. 
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Figure  3.4.  Conventional  and  Capon’s  method  wavenumber  spectral  estimates  for  a  synthetic 
signal  consisting  of  a  sinusoid  with  phase  slowness  of  2  skm'1  propagating  due  south  and  10% 
additive  white  noise. 
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CHAPTER  4 

COMPARISON  OF  AMBIENT  NOISE  WITH 
SURFACE  CONDITIONS 


The  purpose  of  the  CIRCUS  experiment  was  to  gather  noise  data  over  a  wide  fre¬ 
quency  and  wavenumber  band,  and  for  a  long  enough  time  period  to  sample  different  weather 
conditions.  The  one-month  deployment  of  the  array  recorded  significant  variability  in  the  noise 
level.  April  and  May  are  not  times  of  severe  weather  in  the  northern  mid-latitudes  so  there  were 
no  local  storms  during  the  period  of  observation  the  variability  in  the  weather  that  occurred  was 
well  correlated  with  changes  in  the  recorded  noise  levels. 

NOISE  SPECTRA  AND  THE  GSOWM  SWELL  MODEL 

The  characteristics  of  the  noise  level  will  be  described  through  the  power  spectra  of 
the  data  records.  Since  the  maximum  instrument  separation  of  156  m,  the  gross  features  of  the 
spectra  quite  consistent  across  all  the  instruments.  As  a  result,  comparisons  between  noise  lev¬ 
els  and  surface  conditions  will  be  made  with  a  single  instrument,  though  not  necessarily  always 
the  same  one. 

Figures  4. 1-4.4  are  synopses  of  an  entire  OBS  data  set.  Power  spectra  are  computed 
for  each  recording  window  which  occur  at  intervals  of  six  hours.  The  spectra  are  placed  side  by 
side  with  time  increasing  along  the  horizontal  axis,  similar  to  a  sonogram.  Each  vertical  scan 
represent  a  single  recording  window  between  1/32  and  10  Hz  frequency.  The  spectra  are  calcu¬ 
lated  for  data  segments  of  64  seconds  (8192  samples)  by  the  multiple  prolate  spheroidal  func¬ 
tion  taper  method  [Thomson,  1982],  The  data  are  corrected  for  the  instrument  response  to 
acceleration  for  the  case  of  seismometers  or  to  pressure  changes  for  the  hydrophone.  For  the 
purpose  of  plotting,  the  spectra  are  averaged  over  frequency  bands  of  1/12  octaves.  Figure  4.1 
is  data  from  the  vertical  seismometer  of  OBS  8,  Figures  4.2  and  4.3  are  the  two  horizontal  chan- 
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nels  of  OBS  8,  and  Figure  4.4  is  data  from  the  pressure  sensor  of  OBS  6.  The  crystal  hydro¬ 
phone  installed  on  OBS  8  had  poor  response  at  the  lower  frequencies  while  OBS  6  was 
configured  with  a  Cox- Webb  transducer.  However  OBS  6  stopped  recording  after  17  days  days 
on  the  bottom. 

The  are  several  features  of  interest  in  these  plots;  refer  especially  to  Figure  4.1.  The 
main  peak  centered  between  0.1  and  0.2  Hz  visible  at  all  times,  is  the  microseism  peak.  The 
maximum  level  is  variable,  changing  by  a  factor  of  up  to  100  several  times  during  the  month. 
The  frequency  of  the  microseism  peak  is  also  variable  between  0.1  and  0.2  Hz.  At  higher  fre¬ 
quency  a  subsidiary  wide  peak  appears  and  disappears  independently  of  the  level  of  the 
microseism  peak.  There  are  four  time  intervals  —  at  days  116,  120,  129,  and  136  —  where  the 
noise  level  grows  especially  quiet  for  a  period  of  36  hours.  The  horizontal  channels  show  addi¬ 
tional  peaks  at  0.3  and  0.6  Hz,  but  these  are  at  multiples  of  the  microseism  peak  and  may  be 
overtones  in  the  OBS  capsule.  The  variation  in  the  strength  of  the  microseism  peak  in  the  pres¬ 
sure  (Figure  4.4)  is  similar  to  the  vertical  (Figure  4.1)  but  not  identical.  The  variation  at  higher 
frequency  is  not  visible  in  the  pressure  mainly  because  pressure  spectra  have  a  faster  roll-off 
than  ground  acceleration;  pressure  is  more  directly  comparable  to  ground  velocity. 

The  deployment  ships  R/V  Melville  and  US  NS  Narragansett  did  not  remain  on  site 
during  the  recording  period  and  thus  there  are  no  direct  measurements  of  the  surface  environ¬ 
ment.  We  use  instead  the  wind  and  swell  hindcasts  of  the  Global  Spectral  Ocean  Wave  Model 
(GSOWM)  of  the  U.S.  Navy  Fleet  Numerical  Oceanography  Center  (FNOC)  [Clancy  et  al., 
1986].  A  wind-wave  growth  model  is  used  to  compute  the  one-dimensional  swell  spectrum. 
The  energy  is  then  distributed  by  azimuth  according  to  the  fourth  power  of  the  cosine  between 
the  wind  direction  and  the  azimuth.  Propagation  of  the  energy  is  calculated  along  great  circles 
at  the  frequency-dependent  group  velocity.  The  initial  wind  state  is  updated  four  times  a  day 
using  the  most  accurate  wind  data  available.  Wave  energy  is  computed  for  15  frequency  bins 
between  0.039  and  0.308  Hz  and  for  24  points  around  the  compass.  The  directional  spectrum  of 
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the  swell  and  also  the  wind  speed  and  direction  and  some  miscellaneous  parameters  are  output 
for  grid  points  separated  by  2.5°  of  latitude  and  longtitude  at  12  hour  intervals.  Figure  4.5 
shows  the  seven  nearest  grid  points. 

Figure  4.6  is  the  output  of  the  GSOWM  for  32.5°N,  120°W,  the  grid  point  closest  to 
the  CIRCUS  array  site,  at  a  distance  of  60  km.  The  format  of  the  plot  in  Figure  4.6a  is  similar 
to  Figures  4. 1-4.4  except  that  the  spectra  are  average  of  bands  of  1/20  octaves  and  are  at  inter¬ 
vals  of  12  hours. 

A  more  quantitative  comparison  of  the  CIRCUS  noise  data  and  the  GSOWM  will  be 
made  later  but  the  visual  comparison  between  Figure  4.1  and  Figure  4.6a  is  striking.  The  peaks 
in  the  swell  energy  at  days  1 12-115,  121-126,  and  135-142  correspond  to  high  energy  periods  of 
the  microseisms.  Moreover  the  seismic  and  pressure  energy  peak  is  at  twice  the  frequency  of 
the  swell.  The  time  evolution  of  the  swell  peak  frequency  shows  the  deep  water  dispersion 
function.  With  some  imagination  one  can  see  the  evolution  in  the  ocean  bottom  noise  data  as 
well. 

Figure  4.6b  is  a  different  representation  of  the  GSOWM  output.  Each  polar  diagram 
represents  one  time  point.  They  are  offset  horizontally  to  line  up  with  the  corresponding  time  in 
Figure  4.6a.  The  vertical  offset  is  purely  for  clarity  of  display.  The  cumulative  swell  energy  is 
plotted  as  a  function  of  azimuth  with  the  direction  corresponding  to  the  direction  of  wave  pro¬ 
pagation.  The  circles  correspond  to  spectral  energy  of  0.25  and  0.5  m2.  The  wind  direction  and 
speed  arc  also  shown  by  the  arrow.  The  convention  is  that  the  arrow  points  in  the  direction  of 
propagation,  which  is  opposite  to  the  metereological  convention.  The  same  circles  correspond 
to  wind  speeds  of  5  and  10  m  s-1  (about  10  and  20  kt  respectively). 

Another  representation  of  the  wind  speed  and  direction  at  32.5°N,  120°W  appears  in 
Figure  4.7.  The  direction  is  very  steady,  coming  from  NW.  The  speed  varies  between  1  and  1 1 
m  s-1  (2  and  23  kt  respectively).  The  swell  spectral  level  at  0. 1  Hz  and  above  varies  with  the 
wind  speed,  while  the  correlation  between  wind  speed  and  spectral  level  below  0.1  Hz  is  less 
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obvious.  The  wind  speed  and  swell  spectra  at  adjacent  grid  points  differ  in  important  ways 
which  will  be  discussed  further. 

The  accuracy  of  GSOWM  swell  energy  spectrum  can  be  checked  sporadically  with 
wave  height  data  from  the  wave  rider  buoy  at  Begg  Rock  (33°N,  119°W),  at  a  distance  of  110 
km  from  the  nearest  grid  point  and  160  km  from  the  CIRCUS  site.  During  May  1987,  the  data 
link  to  Begg  Rock  was  intermittent  so  the  wave  height  data  are  sparse  (Table  4.1).  The  agree¬ 
ment  between  the  model  and  the  data  varies,  three  examples  are  shown  in  Figure  4.8. 

As  has  been  described  briefly  in  preceding  chapters,  it  has  long  been  recognized  that 
ambient  noise  levels  at  the  ocean  floor  at  frequencies  below  10  Hz  are  primarily  due  to 
processes  at  the  ocean  surface.  The  theory  of  microseisms  will  be  stated  briefly  and  then  a  com¬ 
parison  between  the  ocean  surface  model  and  the  bottom  noise  will  be  made. 


REVIEW  OF  MICROSEISM  THEORY 


The  passage  of  gravity  waves  across  the  water  surface  creates  a  pressure  oscillation 
which  varies  with  depth.  Following  Phillips  [1980],  let  u  be  the  particle  velocity  vector  with 
components  (u,v,w),  and  define  a  potential  such  that  u=V$.  The  boundary  condition  at  the 
ocean  floor  z=-£>  is  d<J>/dz=  0  and  the  surface  displacement  is 


C=a  cos(k-x-co/ ), 

where  k  and  to  are  related  by  the  dispersion  equation 

to2  =  gk  tanh  kD. 

The  velocity  potential  is 


,  c oa  cosh  k(z+D)  .  „  , 

9= - .  .  .  ;n — -sm(k-x-OM). 

xsmh  kD 

The  pressure  to  first  order  as  a  function  of  depth  is  found  from 


,  .  d<}> 

p(z)  =  -pgz-p  -j£- 

Substituting  in  (4.2)  and  (4.3)  yields 


/  \  cosh£(z+Z?)  . 

P  (z )  =  -p  gz  +  pga - —  -~cos(k  -x-cof ). 

cosh  kD 


(4.1) 

(4.2) 

(4.3) 

(4.4) 

(4.5) 
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The  first  term  is  the  hydrostatic  pressure.  The  dynamic  pressure  on  the  bottom  z=-Z),  reduces 
to 


p(D)= 

’  cosh  kD 


(4.6) 


When  Id)  »1  this  quantity  becomes  negligible. 

By  talcing  the  analysis  to  second  order,  it  is  found  that  there  can  be  a  substantial  pres¬ 
sure  signal  on  the  bottom,  even  in  deep  water  [ Miche ,  1940;  Longuet-Higgins,  1950],  Still  fol¬ 
lowing  Phillips  [1980],  the  horizontal  spatial  mean  of  the  pressure  up  to  second  order  tenns  is 


P  (z )  =  gfpdz  +  p(0)  Jy  (0)  -  p(z  )w  (z  )2.  (4.7) 

The  first  term  is  the  hydrostatic  pressure  and  the  last  term  can  be  shown  to  attenuate  exponen¬ 
tially  with  depth.  The  second  term  vanishes  if  the  spatial  mean  is  constant  over  time  as  is  the 
case  for  a  traveling  wave.  However  for  a  standing  wave, 


C  =  a  cos(k-x)cos(cot ), 
this  term  becomes 


(4.8) 


px  =  -yp(0)a2CD2cos2o)f.  (4.9) 

The  frequency  of  the  pressure  oscillation  is  double  that  of  the  surface  wave,  and  the  amplitude 
is  proportional  to  the  square  of  the  amplitude  and  the  square  of  the  frequency  of  the  surface 
wave. 

Longuet-Higgins  [1950]  considered  a  rigid  ocean  bottom,  Hasselman  [1963]  extended 
the  analysis  to  include  an  elastic  medium  below  the  water  layer.  Two  important  results  accrued 
from  this.  First,  the  elastic  medium  was  found  to  be  excited  primarily  by  a  source  in  the  water 
with  phase  velocity  equal  to  that  of  the  resonant  modes  in  the  solid.  Second,  and  as  a  result  of 
the  first,  it  became  necessary  to  consider,  not  standing  waves  with  infinite  phase  velocity  caused 
by  exactly  opposed  wavenumber  components,  but  wave  trains  with  nearly  opposite 
wavenumbers  that  produce  a  pressure  signal  with  finite  phase  velocity.  Hasselman  [1963] 
derived  an  expression  for  the  pressure  spectrum  which  included  terms  with  sums  of 
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wavenumbers  and  differences  of  wavenumbers.  The  sum  terms  result  in  components  with  very 
low  phase  velocity  and  therefore  do  not  contribute  to  the  microseism  spectrum.  The  difference 
terms  become  important  when  kj— k2.  The  approximation  to  the  pressure  spectrum  just  below 
the  ocean  surface  then  becomes 

P (k,co)  =  P2^03  J  Z  (co/2, 0)Z(co/2,0+rc)40.  (4. 10) 

^  —1Z 

The  pressure  oscillation  will  only  propagate  away  from  the  sea  surface  if  Ik  l<co/a.  By  integrat¬ 
ing  over  wavenumber  between  0<Ikl<co/a,  the  frequency  dependence  for  a  point  measurement 
of  the  pressure  on  the  ocean  floor  becomes 

2  2  3  * 

P  (CO)  =  rcp  S  03  f  Z  (co/2 ,0)Z  (co/2, 0+7t)d  0.  (4.11) 

4a2  4 

The  pressure  spectrum  is  determined  by  an  integral  over  azimuth  of  the  product  of  the  direc¬ 
tional  swell  spectrum  with  itself  evaluated  in  the  opposite  direction  and  all  evaluated  at  half  the 
frequency  of  the  pressure.  Equation  (4.11)  includes  no  effect  from  the  bottom. 

Swell  energy  at  a  given  frequency  but  with  opposing  direction  may  have  several 
causes.  Haubrich  et  al.  [1963]  demonstrated  convincingly  that  reflections  of  swell  may  lead  to 
wave-wave  interactions.  There  were  several  lines  of  evidence  for  this  conclusion.  There  was 
no  appreciable  time  delay  between  the  variation  in  the  swell  spectrum  and  the  variation  in  the 
seismic  spectrum  indicating  that  the  generating  area  was  very  close  to  the  coast.  The  measured 
beamwidth  of  the  seismic  energy  was  also  consistent  with  a  generating  area  in  a  narrow  strip  in 
the  ocean  near  the  receiver.  Finally,  the  narrowness  of  the  microseism  spectral  peak  limited  the 
difference  in  propagation  distance  of  the  interacting  swell  sets  and  constrained  the  interaction  to 
be  near  the  coast. 

There  has  not  been  enough  data  collected  to  determine  if  the  microseism  level  in  the 
ocean  varies  with  distance  from  the  coast.  It  is  known  that  there  is  a  double  frequency  microse¬ 
ism  presence  in  the  open  ocean.  An  alternative  way  in  which  wave-wave  interactions  could  be 
generated  is  if  the  wind  direction  varies  rapidly  and  generates  a  new  swell  field  on  top  of  the  old 
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one  [Kibblewhite  and  Ewans,  1986].  However,  though  the  level  may  not  be  as  high,  the  double 
frequency  energy  on  the  ocean  floor  appears  to  exist  independently  of  events  such  as  wind 
shifts.  The  directional  spectra  of  local  wind  waves  appear  to  have  a  spread  of  wavenumbers, 
Longuet-Higgins  et  al.  [1963]  parametrized  the  directional  dependence  of  the  spectrum  of  the 
surface  wave  field  by 


cos^yO-Go)] 

G  (0)  = - - - — • 

2-'-*l7r1/2r(^+i)/r(-y-+|) 


(4.12) 


This  parametrization  was  supported  by  observations  by  Tyler  et  al.  [1974],  They  found  that  the 
beam  parameter  p  varies  between  1  at  high  frequencies  and  10  at  the  spectral  peak. 

If  the  directional  spectrum  of  the  swell  is  not  directly  available,  it  is  often  useful  to 
separate  the  frequency  and  directional  dependence  in  Equation  (4.1 1), 


Z(f,Q)  =  W(f)G(Q)  (4.13) 

[Hughes  1976;  Kibblewhite  and  Ewans,  1986].  In  this  form,  the  pressure  spectrum  becomes 

-  2  4  2  K 

P(f)=  f3W2(f/2)jG  (0)G  (0+*)</  0,  (4.14) 

cc 

where  G  (0)  often  takes  the  form  of  Equation  4.12. 

The  pressure  field  in  turn,  generates  seismic  waves  in  the  ocean  floor,  and  the  passage 
of  seismic  waves  has  an  effect  on  the  pressure.  The  displacement  of  the  horizontal  ocean  floor 
can  be  computed  by  the  convolution  integral 

u(r,z  ,w)  =  J^.G  (r,z  ;r0,z0;co)/  (r-r0^0.“)4  r0,  (4. 15) 

where  /  (r,z,co)  is  a  body  force,  and  G (rrz;r0,z0,co)  is  the  Green’s  function,  the  response  to  a  point 

source  at  (r0,z0).  The  Green’s  functions  can  be  computed  numerically  with  a  variety  of  methods 
[Kennett  and  Kerry,  1979;  Woodhouse,  1980;  Schmidt  and  Tango,  1986].  The  actual  applica¬ 
tion  of  this  formula  is  complicated  in  the  subject  of  microseisms  because  the  source  is  a  par¬ 
tially  correlated  distributed  function  [Hasselman,  1963].  With  some  manipulation,  it  is  possible 
to  recast  Equation  4.15  in  terms  of  a  transfer  function  and  the  spatial  coherence  of  the  source. 
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This  was  used  effectively  by  Adair  [1985]  and  Webb  and  Cox,  [1986].  However  this  will  be  not 
particularly  useful  in  this  study.  The  transfer  function  approach  requires  different  frequencies 
and  wavenumbers  to  be  independent.  The  next  chapter  will  present  data  showing  that  there  is 
scattering  between  wavenumbers  and  that  therefore  the  generation  of  seismic  noise  is  more 
complicated  than  is  allowed  by  the  two-dimensional  Green  function  approach. 

COMPARISON  OF  SWELL  AND  BOTTOM  NOISE 

The  application  of  Equation  (4.11)  to  the  GSOWM  directional  spectra  should  model 
the  ocean  bottom  pressure  spectra.  A  problem  appears  when  this  is  attempted.  In  general,  the 
directional  spectral  generated  by  the  GSOWM  for  this  area  contain  no  energy  directed  in  oppo¬ 
site  directions,  so  application  of  the  integral  over  azimuth  yields  zero.  The  wave  growth  part  of 
the  model  includes  a  cos4  distribution  of  energy  over  azimuth.  However,  the  data  are  provided 
to  three  decimal  digits  and  since  the  quantities  in  the  frequency  and  directional  bins  are  small 
(in  units  of  ft2),  the  information  may  be  lost  in  the  output. 

Instead  of  doing  the  direct  integration  of  Equation  (4.11),  the  one-  dimensional  swell 
spectrum  is  taken  and  multiplied  by  Equation  (4.12)  with  p=2.  Then  Equation  (4.14)  is  applied 
to  generate  a  pressure  spectrum  prediction.  An  example  computed  for  the  grid  point  at  32.5°N, 
120°W  is  shown  in  Figure  4.9.  This  grid  point  provided  a  better  fit  to  the  measurements  than 
the  one  closest  to  the  array  site.  The  latter  point,  though  closer  to  the  measurements  is  at  a  loca¬ 
tion  on  the  continental  borderland  and  is  in  much  shallower  water. 

The  comparison  between  measured  and  predicted  pressure  is  shown  for  three  time 
periods  in  Figure  4.10.  The  degree  of  fit  is  variable.  Certain  high  energy  events  in  the  recorded 
pressure  are  not  predicted  by  the  swell  model.  The  peak  in  the  pressure  spectra  for  days  112- 
1 17  is  much  stronger  than  predicted  and  it  lasts  longer.  As  a  general  rule  there  is  a  local  wind 
wave  component  at  the  higher  frequencies  in  the  GSOWM  which  is  not  visible  in  the  pressure 
spectra.  The  relationship  between  wind  speed  and  swell  spectral  density  at  0.2  Hz  is  clearly 
evident  in  Figure  4.11.  The  wind  state  used  to  compute  wave  growth  in  the  GSOWM  is 
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updated  four  times  a  day  with  actual  measurements  [Clancy  et  al,  1986].  These  grid  points  are 
not  far  from  land  and  major  military  and  commercial  shipping  lanes  so  the  wind  data  is 
expected  to  be  accurate. 

Some  aspects  of  the  OBS  8  vertical  component  fit  the  swell  hindcast  better  than  the 
pressure  does,  especially  at  the  higher  frequencies.  Above  1  Hz  the  seismic  data  exhibit  fairly 
constant  levels  for  a  given  frequency,  interrupted  for  short  intervals  by  troughs  of  very  low 
noise  levels  (Figure  4.1).  The  abrupt  drop  in  noise  level  is  even  more  evident  in  Figure  4.12. 
Examination  of  troughs  shows  that  they  evolve  with  time.  The  falling  edge  occurs  first  at  the 
higher  frequencies  (4  Hz)  and  migrates  down  to  1  Hz  over  the  period  of  about  a  day.  The  rising 
edge  24  to  36  hours  later  follows  the  same  pattern.  These  troughs  occur  at  times  when  the  sur¬ 
face  wind  speed  is  low.  This  step-like  behavior  of  the  spectrum  has  also  been  observed  by 
McCreery  and  Duennebier,  [1989].  When  the  wind  dies  down,  the  waves  at  higher  frequency 
are  attenuated  more  quickly.  The  peak  frequency  of  wind  waves  has  been  observed  to  get  lower 
as  the  duration  of  the  wind  increases  (e.g.  Phillips,  1980),  and  this  is  consistent  with  the  obser¬ 
vations  here.  During  a  sustained  wind  period,  the  wave  spectrum  quickly  reaches  saturation 
with  its  characteristic  /-s  frequency  dependence. 

The  peak  at  days  138-143  is  much  more  prominent  in  the  seismic  data  than  in  the 
GSOWM  at  32.5°N,  120°W.  A  plot  of  a  wider  area  (Figure  4.13)  shows  this  peak  in  the  swell 
more  prominently  to  the  west.  The  predominant  swell  direction  at  this  point  is  from  the  north, 
and  the  GSOWM  propagates  swell  along  great  circles  so  the  grid  points  further  east  may  be  sha¬ 
dowed  by  Pt.  Conception.  Only  one  record  is  available  from  Begg  Rock  at  day  138  and  it  does 
not  differ  radically  from  the  GSOWM  hindcast  (Figure  4.8).  It  appears  then,  that  the  generation 
of  seismic  waves  at  this  time  is  taking  place  far  to  the  west  and  propagating  to  the  array  as  Ray¬ 
leigh  waves.  A  definitive  way  to  test  this  would  be  to  determine  the  direction  of  propagauon  of 
the  seismic  waves  at  this  time.  It  will  be  shown,  however,  in  the  last  chapter  that  the  CIRCUS 
array  was  not  able  to  resolve  velocities  at  the  microseism  peak  frequency. 


Table  4.1 


Bees  Rock  Buov  Records 

no. 

day 

time  (CUT) 

i 

112 

07:05 

2 

112 

19:04 

3 

113 

07:04 

4 

113 

13:04 

5 

114 

07:04 

6 

114 

10:30 

7 

114 

13:04 

8 

114 

19:06 

9 

116 

22:29 

10 

117 

10:29 

11 

118 

07:04 

12 

121 

04:32 

13 

121 

10:30 

14 

121 

13:04 

15 

121 

19:05 

16 

122 

07:04 

17 

123 

19:04 

18 

123 

22:32 

19 

124 

16:30 

20 

124 

19:05 

21 

125 

01:06 

22 

125 

04:26 

23 

125 

07:07 

24 

129 

10:30 

25 

137 

22:30 

26 

139 

01:07 

27 

140 

19:05 
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Figure  4.1.  Spectral  history  of  the  vertical  seismometer  component  of  OBS  8.  The  horizontal 
axis  is  in  days,  starting  with  April  22  and  ending  with  May  23.  Each  vertical  scan  is  a  power 
spectrum  for  a  single  recording  window.  There  are  four  recording  windows  per  day.  The  color 
is  proportional  to  the  decibel  power  as  shown  by  the  color  bar  on  the  right  side.  The  spectra 
have  been  resampled  at  logarithmic  intervals  for  the  purpose  of  plotting. 
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Figure  4.2.  Spectral  history  of  the  first  unrotated  horizontal  component  of  OBS  8.  See  Figure 


4.1  for  details  of  the  plot. 
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Figure  4.3.  Spectral  history  of  the  second  unrotated  horizontal  component  of  OBS  8.  See  Fig¬ 
ure  4. 1  for  details  of  the  plot. 
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Figure  4.4.  Spectral  history  of  the  pressure  component  of  OBS  6.  See  Figure  4.2  for  details  of 
the  plot.  Blank  areas  indicate  missing  data. 
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Figure  4.5.  GSOWM  model  grid  points  near  the  array  site  for  which  data  was  made  available 
by  the  Fleet  Numerical  Oceanography  Center. 
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Figure  4.6.  Spectral  history  of  the  GSOWM  swell  model,  a)  Spectral  power  as  a  function  of  fre¬ 
quency,  summed  over  azimuth  and  displayed  the  same  as  in  Figure  4.1,  b)  Directional  power 
summed  over  all  frequencies.  The  circles  mark  spectral  energy  levels  of  0.25  and  0.5  m2.  The 
origin  of  the  circles  for  each  time  coincides  with  the  corresponding  time  in  4.6a.  The  wind  vec¬ 
tor  is  also  drawn  for  each  time. 
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Figure  4.7.  GSOWM  wind  speed  and  direction  hindcast  for  the  nearest  grid  point  to  the 
CIRCUS  array. 
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Figure  4.8.  Comparison  of  the  GSOWM  swell  hindcast  spectrum  with  data  from  the  buoy  at 
Begg  Rock  at  three  times.  The  solid  lines  are  for  the  Begg  Rock  data  computed  from  one- 
minute  records.  The  dotted  lines  are  the  swell  model,  a)  Day  1 13,  b)  Day  1 17,  c)  Day  138. 
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Figure  4.9.  Comparison  of  the  bottom  pressure  spectrum  predicted  from  the  GSOWM  swell 
model  and  the  CIRCUS  measurements.  The  form  of  the  plot  is  similar  to  those  of  Figures  4.6a. 


a)  GSOWM  hindcast,  b)  CIRCUS  data. 
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Figure  4.10.  Comparison  of  the  bottom  pressure  spectrum  predicted  from  the  GSOWM  swell 
model  and  the  CIRCUS  measurements  at  three  times.  The  solid  lines  are  the  CIRCUS  measure¬ 


ments.  The  dotted  lines  are  the  prediction  from  the  GSOWM  swell  hindcast.  a)  Day  113,  b) 
Day  120,  c)  Day  126. 
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CSOWM  Wind  Speed 


Figure  4.1 1.  Comparison  of  GSOWM  hindcasts  for  wind  and  swell  energy  at  0.2  Hz.  a)  Wind 
speed,  b)  Swell  energy  at  0.2  Hz. 
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FNOC  Wind  Hindcast  at  32.5°N.120°W 


Figure  4.12.  Comparison  of  the  GSOWM  hindcast  for  wind  speed  and  the  measured  seismic 
noise  energy  at  2  Hz.  a)  Wind  speed,  b)  Seismic  level  at  2  Hz. 
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Figure  4.13.  GSOWM  swell  hindcasts  at  several  grid  points  near  the  array.  See  Figure  4.6a  for 


description  of  the  plot. 
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CHAPTER  5 

COHERENCE  MEASUREMENT  AND  MODELING 


Multiple  sensors  allow  the  direct  measurement  of  signal  coherence,  and  if  sensors  are 
spaced  at  a  variety  of  offsets,  coherence  length  can  be  estimated.  The  coherence  length  of 
ambient  noise  has  two  important  consequences  for  arrays  designed  to  sense  known  signals.  The 
first  is  that  the  array  gain  only  becomes  substantial  if  the  sensor  separations  are  greater  than  the 
noise  coherence  length.  The  second  is  that  differences  in  coherence  lengths  (or  wavelengths) 
can  be  used  to  improve  the  signal  to  noise  ratio.  Here  the  noise  has  short  wavelengths  and  sen¬ 
sor  groups  act  as  wavenumber  filters  to  suppress  short-wavelength  noise  while  leaving  longer- 
wavelength  signals  (from  ships,  earthquakes,  or  explosions)  unimpaired.  This  technique  is 
widely  used  in  sensor  arrays  used  in  reflection  seismology  both  on  land  and  at  sea. 

The  nine  instruments  give  36  different  sensor  separations.  All  the  separations  greater 
than  80  m  are  relative  to  OBS  1  which  was  separated  from  the  rest  of  the  array.  The  coherence 
for  each  pair  of  instruments  for  a  particular  window  are  shown  in  Figure  5.1.  The  plot  is  typical 
of  the  observations  at  this  site.  There  are  two  bands  of  high  coherence,  below  0.4  Hz  and 
between  0.8  and  3  Hz.  Between  those  two  bands  there  is  a  notch  of  lower  coherence  which 
widens  with  increasing  separation. 

The  ambient  seismic  noise  field  derives  its  characteristics  from  three  processes:  the 
source  of  energy,  the  coupling  of  the  energy  into  the  elastic  medium,  and  the  propagation  of 
energy  in  the  elastic  medium.  It  was  shown  in  Chapter  4,  through  the  theory  of  Longuet- 
Higgins,  [1950]  and  Hasselman  [1963],  that  under  certain  conditions,  ocean  surface  swell  can 
generate  pressure  fluctuations  that  reach  the  ocean  floor  unattenuated.  It  was  also  shown  that 
the  wavenumber  of  the  resulting  pressure  fluctuations  is  limited  to  I  k  I  <w/a,  where  a  is  the 
acoustic  velocity  of  water.  The  pressure  fluctuations  can  couple  into  the  ocean  floor  when  their 
wavenumber  is  equal  to  that  of  the  resonant  modes  in  the  elastic  medium.  However,  the 
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excitation  of  high  wavenumber  (low  velocity)  seismic  waves  which  are  observed  requires  a  dif¬ 
ferent  interaction  which  is  not  well  understood.  The  last  process,  the  propagation  of  energy  in 
the  elastic  medium,  will  be  the  main  focus  of  this  chapter.  The  observed  spatial  coherence  puts 
strong  constraints  on  the  nature  of  the  seismic  propagation  between  the  coupling  region  and  the 
receivers. 

The  coherence  model  is  calculated  from  seismic  modes  in  the  sedimentary  wave 
guide.  Modal  dispersion  curves  and  synthetic  seismograms  were  calculated  using  the  normal 
mode  formalism  of  Woodhouse  [1980]  as  implemented  by  Gomberg  and  Masters  [1988].  The 
mode  spectra  are  calculated  for  a  point  tensor  source  in  the  lossy,  elastic  solid. 

Measured  seismic  velocities  are  used  in  the  input  model.  A  1984  experiment,  under¬ 
taken  at  DSDP  hole  469,  studied  the  velocity  structure  of  the  sediments  [Sauter  et  al.,  1986]. 
Explosives  were  detonated  on  the  ocean  floor  so  that  significant  energy  would  be  coupled  into 
the  sediment  waveguide.  Dispersed  Stoneley  wave  trains  were  generated  and  the  observed 
dispersion  pattern  (Figure  5.2)  was  matched  by  iterative  modeling  with  the  synthetic  seismo¬ 
gram  program.  The  shape  of  the  envelope  of  the  wavetrains  is  strongly  determined  by  the 
seismic  Q  structure  of  the  waveguide.  A  Q  model  was  determined  by  matching  the  synthetic 
seismograms  to  the  data  [Sauter,  1987].  The  bulk  of  the  energy  of  these  dispersed  wave  trains 
propagated  at  velocities  between  30  and  100  m  s-1  and  was  quickly  attenuated. 

The  observation  that  spectral  energy  from  small  explosions  in  shallow  water  produced 
Stoneley  waves  with  a  spectral  peak  in  the  1-5  Hz  range,  led  Tuthill  et  al.,  [1981]  to  suggest  that 
these  waves  were  significant  contributors  to  sea-floor  noise.  The  CIRCUS  experiment  was 
designed  to  examine  this  assertion  by  attempting  to  observe  natural  occurrences  of  these  very 
low  velocity  waves  on  the  ocean  floor. 

The  present  study  is  concerned  with  frequencies  extending  much  lower  than  were 
observed  from  the  small  explosions.  The  wavefunction  at  the  lower  frequencies  extends  into 
the  basaltic  crust  so  the  sediment  velocity  model  of  Sauter  et  al.  [1986]  has  been  merged  with  a 
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crustal  model  which  is  representative  of  this  region  of  the  Pacific  ocean  crust  [Spudich  and 
Orcutt,  1980].  The  velocity  model  is  shown  in  Figure  5.3  and  the  associated  phase  and  group 
dispersion  curves  in  Figure  5.4.  The  dispersion  model  is  plotted  as  a  function  of  slowness 
because  this  parameter  relates  more  directly  to  coherence  length  than  does  phase  velocity.  The 
phase  velocity  undergoes  a  change  at  0.5  Hz  as  it  makes  the  transition  between  Rayleigh  waves 
and  Stoneley  waves.  The  Stoneley  (or  Scholte)  waves  are  the  continuation  of  Rayleigh  waves 
to  phase  velocities  lower  than  that  of  water,  1.5  km  s"1  (c.f.  Tolstoy,  1954).  At  these  lower  velo¬ 
cities,  the  vertical  wavefunction  is  evanescent  in  the  water,  decaying  upward  away  from  the 
interface.  The  sea  floor  becomes  the  surface  around  which  the  waves  are  trapped. 

Computation  of  plane  wave  synthetic  seismograms  alone  makes  no  predictions  for  the 
coherence  structure,  however  Capon  [1973]  showed  that  limits  on  coherence  length  could  be 
imposed  by  an  isotropic  distribution  of  sources.  For  an  isotropic  distribution  of  plane  wave 
sources,  Capon  [1973]  showed  by  integrating  the  wavefield  over  azimuth  that  the  coherence 
field  is  given  by 

tf2(m,0  =  7o(<op((o)r)  (5.1) 

where  K2  is  the  mean  square  coherence  (MSC),  J0  is  the  zero-order  Bessel  function  of  the  first 
kind,  to  is  the  radial  frequency,  p(co)  is  the  frequency  dependent  phase  slowness  (the  reciprocal 
of  phase  velocity),  and  r  is  the  separation  between  sensors.  The  function  p( to)  in  the  argument 
to  the  Bessel  function  is  determined  from  the  phase  dispersion  functions  (Figure  5.4)  A  realiza¬ 
tion  of  Equation  (5.1)  for  mode  0  is  shown  in  Figure  5.5.  The  distance  ordinate  of  the  first  zero 
of  the  Bessel  function,  by  which  coherence  length  is  being  defined,  is  directly  related  to  the 
phase  slowness  at  a  given  frequency.  The  transition  from  long  to  short  coherence  length  occurs 
between  0.5  and  1  Hz  and  is  associated  with  phase  velocity  changes  in  this  frequency  band. 

In  contrast  to  Figure  5.5,  the  observations  (Figure  5.1)  display  multiple  transitions 
from  long  to  short  coherence  length,  evidence  that  multiple  modes  are  present.  The  mean 
squared  coherence  measures  the  proportion  of  energy  that  is  organized.  The  total  squared 
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coherence  is  the  sum  of  the  squared  coherence  from  many  sources  in  proportion  to  the  power  of 
each  individual  source.  For  the  present  model, 

K\ CO)  =  2>;(co)  Jo  (cop,(co)r)  (5.2) 

where  wj  is  the  proportional  power  in  mode  j. 


w,(co)  = 


Pjm 


i-O 


and  Pj(n)  is  the  power  in  mode  j. 

Power  spectra  for  modes  0  through  3  are  shown  in  Figure  5.6,  computed  for  a  distance 
of  15  km  and  a  source  depth  of  1  m  below  the  sea  floor.  The  higher  modes  dominate  at  frequen¬ 
cies  above  1  Hz.  Figure  5.7  shows  the  coherence  predicted  from  Equation  (5.2)  using  the  pro¬ 
portional  power  from  Figure  5.6.  Comparison  of  the  coherence  model  (Figure  5.7)  and  the 
coherence  observations  (Figure  5.1)  reveal  a  good  match.  The  computation  explains  the  band 
of  high  coherence  between  0.8  and  3  Hz  and  the  low  coherence  notch  between  0.4  and  0.8  Hz. 

To  compute  synthetic  power  spectra,  a  source  distance  and  depth  must  be  assumed.  A 
distance  of  15  km  was  chosen  by  trial  and  error  to  achieve  the  observed  separation  between  the 
bands  of  high  coherence.  The  the  lower  edge  of  the  high  frequency  band  in  the  coherence 
model  (Figure  5.7)  is  controlled  by  the  frequency  at  which  the  higher  modes  begin  to  dominate 
the  fundamental  mode.  The  eigenfunction  for  the  mode  0  is  strongest  near  the  ocean  floor. 
Mode  0  is  thus  excited  strongly  by  a  source  near  the  surface,  but  since  the  Q  of  the  upper  sedi¬ 
ments  is  low,  mode  0  is  attenuated  more  strongly  than  are  the  higher  modes.  This  is  illustrated 
by  four  examples  in  Figure  5.8  for  different  ranges  and  source  depths.  For  a  shallow  (1  m) 
source  the  dependence  of  mode  attenuation  on  distance  is  strong.  For  a  deeper  (100  m)  source 
the  mode  attenuation  with  distance  is  less  clear.  In  Figure  5.9  the  frequency  at  which  the  higher 
modes  exceed  the  fundamental  mode  is  plotted  as  a  function  of  source  distance  and  depth.  For 
a  deeper  source,  the  fundamental  mode  is  not  excited  as  strongly  and  there  is  less  dependence 
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a  deeper  source,  the  fundamental  mode  is  not  excited  as  strongly  and  there  is  less  dependence 
on  distance.  The  source  moment  tensor  may  also  affect  the  excitation  of  modes.  Since  we  have 
insufficient  data  to  determine  such  a  tensor  a  diagonal  tensor  has  been  used  for  these  calcula¬ 
tion. 

A  physical  interpretation  can  be  associated  with  the  15  km  distance  derived  from  the 
coherence  model.  Figure  5.10  shows  a  bathymetric  map  for  the  immediate  vicinity  of  the 
experiment,  with  a  circle  marking  a  radius  of  15  km  from  the  array  site.  Significant  topographic 
features  are  located  roughly  at  this  distance,  in  particular  the  Patton  escarpment  to  the  east  and  a 
seamount  to  the  northwest  A  likely  explanation  is  that  inhomogeneities  in  the  boundaries  can 
scatter  incident  waves  into  the  surficial  layers.  Levander  and  Hill  [1985]  showed  by  numerical 
modeling  how  high  wave  number  signals  generated  at  depth  (or  analogously  at  the  water  sur¬ 
face)  can  be  coupled  into  high  wave  number  (slowly  propagating)  signals  in  the  ocean  bottom. 
Dougherty  and  Stephen  [1988]  demonstrated  how  volume  heterogeneity  could  convert  incident 
body  waves  into  surface  waves.  While  we  cannot  make  a  definitive  distinction  here  between 
these  two  possibilities,  the  bathymetric  data  available  to  us  most  closely  resembles  the  assump¬ 
tions  made  by  Levander  and  Hill  [1985]. 

Low  coherence  near  0.6  Hz  in  this  model  requires  that  the  coherence  between  the 
separate  modes  is  small.  All  modes  are  excited  by  the  same  source  but  because  of  the  large 
variation  in  phase  velocity,  different  frequencies  reaching  the  receiver  at  a  given  time  were  gen¬ 
erated  at  different  times.  As  an  extreme  example,  mode  3  at  1  Hz  will  propagate  15  km  in  10  s, 
while  mode  0  at  1  Hz  requires  150  s.  The  time  difference  is  long  compared  to  the  correlation 
time  of  ocean  swell.  The  expected  time  separation  between  modes  exceeds  the  length  of  the 
discrete  Fourier  transform  windows  and  the  calculated  coherence  is  therefore  reduced  (c.f.  Jen¬ 
kins  and  Watts,  1968,  p.  399). 

The  coherence  model  differs  from  the  observed  coherence  in  the  following  respects. 
The  low  coherence  notch  is  not  nearly  as  deep  in  the  data  and  the  band  of  coherence  near  1  Hz 
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is  wider  in  the  data,  which  suggests  that  the  coherence  lengths  for  the  model  are  too  low.  A 
more  general  beam  pattern  than  for  the  isotropic  case,  is  the  form  cos'1  (0/2)  [Longuet-Higgins  et 
al.,  1963;  Webb,  1986].  An  isotropic  wave  field  is  implied  by  p  =  0,  and  a  single  plane  wave  by 
p  =  oo.  Coherence  functions  are  shown  in  Figure  5.11  for  a  variety  of  beam  parameters  p.  The 
coherence  length  increases  for  increasing  directional  confinement,  suggesting  that  a  narrower 
beam  can  improve  the  fit.  The  coherence  model  of  Figure  5.7  is  recalculated  for  p  =  4  and 
shown  in  Figure  5.12  and  shows  a  better  match  to  the  data.  The  improvement  with  a  narrower 
source  beam  is  not  surprising  since  the  noise  is  probably  associated  with  the  dominant  bathy¬ 
metric  feature,  which  is  the  Patton  escarpment  to  the  northeast. 

The  dimension  of  the  array  is  very  small  compared  to  the  distance  to  the  source.  If 
the  array  had  been  much  more  extensive,  so  that  we  could  decompose  the  wavefield  by  mode 
and  azimuth,  we  could  extract  more  information  about  the  scattering  sources  and  mechanisms. 
In  the  study  of  earthquakes,  the  global  array  of  low  frequency  seismic  stations  allows  recon¬ 
struction  of  the  time  history  of  each  component  of  the  tensor  source  [Dziewonski  et  al.,  1987; 
Riedesel  and  Jordan,  1989;  Gilbert,  1971],  In  principle  the  only  difficulty  barring  us  from 
doing  this  at  smaller  scale  on  the  sea  floor  is  the  matter  of  achieving  adequate  spatial  sampling 
to  decompose  the  wavefield. 

In  summary,  measurement  with  a  150  m  aperture  array  on  the  seafloor  has  placed 
significant  constraints  on  the  coherence  of  ambient  noise  between  0.05  and  10  Hz.  The 
observed  pattern  of  coherence  is  well  modeled  by  an  isotropic  or  broad  beam  distribution  of 
low  order  Raleigh/Stoneley  waves.  The  velocity  and  attenuation  structure  of  the  Rayleigh  wave 
guide  affects  the  coherence  structure  by  controlling  the  relative  power  of  the  different  modes  as 
a  function  of  frequency  and  distance  from  the  source.  The  coherence  field  at  this  site  is  pro¬ 
foundly  affected  by  the  low  velocity  sediments  that  cover  the  basin.  In  the  framework  of  this 
model,  the  band  of  high  coherence  length  between  0.8  and  2  Hz  suggests  that  the  seismic  waves 
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which  comprise  the  noise  are  excited  at  a  distance  of  15  to  20  km  or  alternatively  at  a  depth 
greater  than  100  m  below  the  ocean  floor.  Either  of  these  possibilities  imply  that  the  energy 
which  is  generated  at  the  ocean  surface  is  scattered  into  the  sedimentary  layer  at  the  boundaries 
of  the  basin.  The  good  fit  of  this  model  to  the  observed  coherence  indicates  that  at  this  site  at 
least,  scattered  energy  is  a  more  important  source  of  seismic  noise  at  the  ocean  floor  than  is 
direct  radiation  from  the  sea  surface  above  the  sensor. 

In  contrast  with  the  results  presented  here,  the  physics  by  which  noise  is  coupled  into 
the  sea  floor  in  shallow  water  is  quite  different  There  Schmidt  and  Kuperman  [1988]  have 
demonstrated  that  it  occurs  by  direct  coupling  of  sea-surface  noise  into  the  waveguide  through 
the  evanescent  part  of  the  wave  function.  Energy  in  the  water  column  has  a  much  stronger  con¬ 
tribution  to  noise  levels  measured  on  the  deep  ocean  floor. 
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Figure  5.1.  Spatial  coherence  measured  from  record  229.  The  horizontal  axis  is  frequency  from 
0.03  to  10  Hz,  and  the  vertical  axis  is  separation  between  sensors  from  0  to  156  m.  Color  is  pro¬ 
portional  to  mean  squared  coherence. 
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Figure  5.2.  Stoneley  wave  recorded  from  an  ocean  floor  shot  The  dispersion  analysis  below  the 
time  series  shows  the  variation  in  frequency  with  time  relative  to  the  source,  or  equivalently,  to 
group  slowness. 
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Figure  5.3.  Velocity  model  for  the  CIRCUS  array  site.  The  velocities  in  the  upper  sediments 
are  derived  from  the  dispersion  analysis  of  bottom  shots.  The  crustal  velocities  are  from  Spud - 
ich  and  Orcutt  [1980]. 
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Dispersion  Model 


Figure  5.4.  Synthetic  dispersion  curves.  The  line  style  is  consistent  for  each  mode  in  both 


displays,  a)  Group  slowness,  b)  Phase  slowness. 
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Figure  5.5.  Synthetic  spatial  coherence  calculated  from  the  mode  0  dispersion  curve.  See  Fig¬ 


ure  5.1  for  description  of  the  plot. 
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Figure  5.6.  Synthetic  power  spectra  and  relative  power  for  modes  0  through  3.  The  line  styles 
are  the  same  as  in  Figure  5.4.  a)  Power  spectra,  b)  Proportion  of  each  mode  to  total  power. 
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Figure  5.7.  Synthetic  spatial  coherence  for  modes  0  through  3  and  an  isotropic  source  distribu¬ 
tion. 
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Figure  5.8.  Examples  of  mode  spectra  for  different  source  depth  and  distances,  a)  Source  depth 
=  1  m,  distance  =  10  km,  b)  Source  depth  =  1  m,  distance  =  100  km,  c)  Source  depth  =  100  m, 
distance  =  10  km,  c)  Source  depth  =  100  m,  distance  =  100  km. 
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Cross- over  frequency  vs.  distance 


Figure  5.9.  Frequency  at  which  one  of  the  higher  modes  becomes  dominant  over  the  fundamen¬ 
tal  mode  versus  source  distance  at  source  depth  of  10  and  100  m. 


Figure  5.10.  Bathymetry  of  the  CIRCUS  array  vicinity.  The  contour  interval  is  100  m. 
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Coherence  Length  vs.  Beam  Parameter 


Figure  5.1 1.  Coherence  versus  distance  for  four  beam  parameters. 
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Figure  5.12.  Synthetic  spatial  coherence  for  modes  0  through  3  and  a  distributed  source 
predominantly  on  one  side  of  the  array. 
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CHAPTER  6 

BEAMFORMING  RESULTS 


Sensors  are  placed  in  arrays  for  two  reasons,  to  enhance  signal  to  noise  ratio,  and  to 
make  direct  measurement  of  phase  velocity,  which  can  be  converted  to  bearing  and  velocity 
estimates.  Arrays  have  been  used  very  successfully  in  the  disciplines  of  radio  science  and 
underwater  acoustics  to  determine  the  bearing  of  signal  sources.  In  both  those  fields,  the  prob¬ 
lem  is  simplified  because  acoustic  waves  in  water  and  electromagnetic  waves  propagate  at 
approximately  a  single,  well  known  velocity.  Apparent  phase  velocity  across  the  array  is  thus 
mapped  uniquely  into  angle  of  incidence.  Seismic  waves,  on  the  other  hand,  travel  at  a  range  of 
velocities  between  40  m  s-1  in  water-saturated  muds  on  the  ocean  floor,  to  13  km  s-1  in  the  lower 
mantle.  Some  additional  knowledge  of  the  propagation  path  is  usually  required  in  the  seismic 
case.  Exploration  of  Earth  structure,  of  the  type  done  for  petroleum  exploration  routinely  uses 
linear  arrays  of  many  elements,  both  for  estimating  apparent  phase  velocity  through  moveout 
determination,  and  for  enhancing  coherent  arrivals  through  stacking.  The  technique  has  not 
been  used  as  commonly  for  global  seismology  and  for  natural  seismic  sources  because  of  the 
more  expensive  sensors  required  for  those  applications.  There  are  a  few  examples  of  large 
arrays  used  for  measuring  teleseismic  signals,  the  most  notable  example  is  LASA,  the  Large 
Aperture  Seismic  Array  [ Capon  etal.,  1969]. 

Data  from  LASA  had  a  substantial  impact  on  the  study  of  microseisms  [Haubrich  and 
McCamy,  1969;  LaCoss  et  al.,  1969].  It  was  determined  that  background  noise  propagates  at 
3.5  km  s'1  in  the  frequency  range  of  0.03  and  0.05  Hz.  At  these  frequencies,  this  velocity  can  be 
identified  with  the  fundamental  Rayleigh  mode  for  a  continental  path.  The  apparent  direction  to 
the  source  of  this  energy  was  from  the  coast  nearest  to  oceanic  storms,  thus  corroborating  the 
hypothesis  that  ocean  waves  are  the  primary  source  of  microseisms.  In  the  frequency  range  of 
0.3  to  0.5  Hz,  the  propagation  was  also  found  to  be  3.5  km  s-1,  corresponding  to  higher  mode 


103 


104 


Rayleigh  waves.  At  frequencies  above  3  Hz,  the  velocity  was  measured  at  13  km  s'1  and  was 
associated  with  P-waves  propagating  in  the  lower  mantle. 

Few  attempts  at  beamforming  on  the  deep  ocean  floor  have  been  made.  The  CIRCUS 
was  the  first  that  permitted  coherent  processing  of  signals  on  more  than  two  sensors.  The 
CIRCUS  array  was  of  a  more  limited  nature  than  typical  land  arrays.  Ocean  bottom  seismome¬ 
ters  are  expensive  and  few  in  number  so  it  was  necessary  to  limit  the  wavenumber  bandwidth  of 
the  array  to  between  6  and  125  cycles  km-1.  The  width  of  the  central  lobe  of  the  array  response 
is  10  km-1,  or  \0l2nf  skm-1  in  terms  of  slowness  (reciprocal  velocity).  Figure  6.1  shows  how 
this  resolution  limit  relates  to  expected  phase  velocities  from  the  dispersion  curves.  Measure¬ 
ments  below  the  the  dashed  line  cannot  be  resolved  with  conventional  beamforming  from  a 
slowness  of  zero.  Some  improvement  in  resolution  can  be  achieved  with  Capon’s  method  but 
the  array  can  effectively  measure  phase  velocity  only  above  0.5  Hz.  The  uncertainty  in  relative 
timing  of  the  OBS’,  which  was  discussed  in  Chapter  2,  will  further  limit  the  resolution  of  beam¬ 
forming. 

Another  limitation  is  that  the  CIRCUS  array  was  never  effectively  calibrated.  The 
experiment  plans  called  for  a  series  of  shots  to  be  detonated  in  the  water  at  a  variety  of  azimuths 
around  the  array.  However  the  OBS  clocks  and  the  clock  by  which  the  shooting  was  timed 
were  offset  from  each  other  by  one  minute,  so  the  first  arrivals  from  the  shots  were  not  recorded 
by  the  OBS.  The  deployment  ships  provide  a  less  satisfactory  alternative  signal  to  be  used  for 
array  calibration.  Rotating  machinery  in  the  propulsion  systems  generates  spectral  lines  in  the 
water  (Figure  6.2).  Impulsive  sources  are  more  desirable  because  timing  and  position  errors  can 
be  seen  directly  on  each  sensors’  record  as  time  offsets.  With  continuous  wave  sources  such 
errors  are  only  observable  as  phase  delay,  and  there  is  an  inherent  ambiguity,  which  becomes 
harder  to  overcome  at  the  higher  frequencies  associated  with  ship  noise.  There  were  two  ships 
present  during  the  first  two  days  of  the  deployment,  the  R/V  Melville  and  the  USNS  Narragan- 
sett.  They  have  radically  different  propulsion  systems  and  so  would  be  expected  to  have  distin- 
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guishable  acoustic  signatures.  It  has  however,  proved  difficult  to  assign  spectral  lines  to  one 
ship  or  the  other.  There  is  usually  no  simple  relationship  between  between  spectral  line  ampli¬ 
tude  and  ship  distance.  Beamforming  has  provide  unambiguous  answers  for  only  two  spectral 
lines,  examples  of  which  are  shown  in  Figures  6.3  and  6.4.  In  these  figures,  the  wavenumber 
axes  have  been  scaled  by  the  analysis  frequency  so  that  the  units  become  those  of  slowness. 
The  vector  fro  the  origin  to  any  point  represents  propagation  in  that  direction  with  slowness 
proportional  to  the  length  of  the  vector.  The  dashed  circle  centered  at  the  origin  represents  the 
acoustic  slowness  of  water  1/a.  A  signal  propagating  horizontally  horizontally  in  the  water  will 
be  indicated  by  wavenumber  analysis  as  a  point  on  this  circle.  A  signal  propagating  vertically 
will  appear  at  the  center.  The  color  at  each  point  represents  the  power  level.  The  5.75  Hz  line 
has  been  identified  with  the  R/V  Melville  while  it  is  holding  station  nearly  above  the  array  (Fig¬ 
ure  6.3).  The  line  is  absent  when  the  ship  is  underway  or  at  some  distance.  The  19.5  and  19.6 
Hz  lines  have  been  identified  with  a  continuous  wave  source  that  was  towed  behind  the  USNS 
Narragansett  (Figure  6.4).  However,  in  the  figure  it  is  evident  that  at  this  frequency  there  is 
also  a  source  with  the  bearing  corresponding  to  the  R/V  Melville.  Figure  6.4  show  a  misfit  of 
the  predicted  position  of  0.2  s  km-1.  In  many  other  cases,  beamforming  does  not  identify  a  bear¬ 
ing  that  corresponds  to  either  one  of  the  ship  positions.  It  is  unknown  at  this  time  whether  the 
inconsistent  performance  of  the  wavenumber  processing  is  due  to  errors  in  the  data  such  as  tim¬ 
ing  errors  or  position  errors,  or  whether  it  is  due  to  low  stability  of  the  spectral  estimates. 

Keeping  in  mind  the  possible  inaccuracies  in  the  wavenumber  estimates,  a  survey  of 
the  early  part  of  the  experiment  follows.  The  error  in  the  clocks  is  expected  to  get  worse  with 
time,  and  indeed  more  scatter  in  the  estimated  velocities  is  evident  in  the  later  records.  Figures 
6.5  through  6.7  show  wavenumber  spectra  for  window  219  for  frequencies  of  0.25,  0.75,  and 
1.25  Hz.  The  spectra  are  computed  by  Capon’s  method  with  segment  averaged  Fourier 
transforms  with  a  time -bandwidth  product  of  25.  Energy  propagating  to  the  southwest  at  each 
frequency.  A  slowly  propagating  wave  is  detected  at  0.75  Hz,  but  otherwise  the  observed  slow¬ 
ness  is  between  0.5  and  1.0  skrn"1.  As  expected,  examination  of  other  records  shows  some 
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variability,  but  slow  propagation  velocity  is  commonly  observed  at  frequencies  between  0.25 
and  0.75  Hz.  Figure  6.8  is  a  summary  of  all  velocities  measured  for  the  first  20  records  between 
0.125  and  3  Hz.  Also  shown  in  the  figure  are  the  dispersion  curves  from  Sauter  et  al.  [1986] 
which  were  used  in  Chapter  5  to  compute  synthetic  spectra.  The  dispersion  curves  are 
emphasized  where  the  individual  mode  is  dominant  (compare  with  Figure  5.4).  There  is  rough 
agreement  with  the  expected  velocities  so  the  results  from  the  beamforming  analysis  are  con¬ 
sistent  with  the  results  derived  from  modeling  the  spatial  coherence.  The  distribution  of  energy 
among  the  modes  is  evident.  As  frequency  increases,  the  phase  slowness  increases  as  it  follows 
a  given  mode’s  dispersion  curve,  then  decreases  when  it  jumps  to  the  next  higher  mode  as  it 
becomes  predominant.  The  scatter  of  the  points  is  relatively  large.  However,  comparison  with 
the  results  of  Haubrich  and  McCamy  [1969]  (see  their  Figure  9),  reveals  that  even  at  LASA, 
under  more  ideal  circumstances,  the  scatter  is  substantial. 

The  direction  of  propagation  is  summarized  in  Figure  6.9.  The  histogram  shows  a 
preponderance  of  energy  propagating  to  the  southwest  This  is  consistent  with  a  source  of 
energy  at  the  edge  of  the  sedimentary  basin,  which  is  closest  to  the  array  site  towards  the 
northeast,  at  the  base  of  the  Patton  escarpment. 

Qear  from  Figures  6.8  and  6.1  is  that  the  target  wavenumber  band  (and  thus  slow¬ 
ness)  of  the  CIRCUS  array  was  too  low.  It  had  been  presumed  from  observation  of  the  close 
range  bottom  shots,  that  the  predominant  energy  would  be  on  the  very  slow  part  of  the  disper¬ 
sion  curve,  above  10  skm-1.  The  discovery  that  the  energy  is  coupled  into  the  sediment  wave 
guide  at  a  distance  of  several  kilometers  means  that  these  very  slow  waves  are  attenuated  below 
detection  by  the  time  they  reach  the  array.  At  the  frequency  at  which  the  fundamental  mode 
would  be  propagating  at  10  to  30  s  km-1,  the  faster  overtones  are  detected  instead. 

The  most  important  result  from  this  study  is  that  there  is  significant  energy  propagat¬ 
ing  at  the  ocean  bottom  at  low  velocity,  especially  at  frequencies  between  0.5  and  3  Hz.  The 
source  of  ambient  noise  energy  is  pressure  fluctuations  at  the  ocean  surface.  It  has  been  shown 
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that  this  type  of  source  cannot  excite  modes  in  the  ocean  bottom  at  phase  slowness  greater  than 
l/a  =  0.66  s  km-1.  The  fact  that  such  low  velocities  are  indeed  observed  indicates  that  some  form 
of  energy  scattering  is  taking  place.  The  nature  of  the  scattering  is  beyond  the  scope  of  this 
study,  but  a  reasonable  conclusion  is  that  it  occurs  where  the  sediment  bed  pinches  out  and 
there  is  a  sloping  boundary  between  the  sediments  and  the  hard  rock. 

IMPLICATIONS  FOR  FUTURE  WORK 

The  decomposition  of  the  ambient  noise  field  into  seismic  modes  has  important  impli¬ 
cations  for  sub-bottom  sensor  data.  The  justification  given  for  making  the  effort  to  deploy  sen¬ 
sors  in  a  borehole  is  that  there  is  a  potential  for  improving  the  signal  to  noise  ratio  of  measured 
teleseismic  signals.  The  degree  to  which  this  is  true  is  a  consequence  of  the  depth  dependence 
of  the  ambient  noise  compared  to  that  of  the  desired  signal.  Body  waves  are  often  of  interest, 
especially  at  higher  frequencies,  and  they  do  not  depend  on  distance  from  the  seafloor  directly. 
Ambient  noise  is  presumed  to  consist  of  interface  waves,  an  assertion  supported  by  this  study. 
The  depth  dependence  of  interface  waves  can  be  determined  from  the  mode  structure.  At  this 
site,  the  fundamental  Rayleigh  mode  is  the  primary  contributor  to  noise  below  0.8  Hz.  Above 
that,  and  up  to  about  5  Hz,  higher  Rayleigh  modes  have  higher  amplitudes.  The  higher  modes 
decay  less  rapidly  with  depth  than  the  fundamental  so  it  is  important  to  consider  their  effects 
when  considering  the  question  of  sensor  depth. 

An  alternative  to  placing  sensors  in  a  borehole,  is  to  use  an  array  of  sensors  to 
enhance  the  signal  to  noise  ratio.  Knowledge  of  the  coherence  length  of  the  ambient  noise  at 
different  frequencies  makes  possible  the  design  of  optimum  array  configurations  to  maximize 
array  gain.  The  present  study  showed  that  the  minimum  coherence  length  in  the  frequencies 
observed  is  50  m,  and  although  not  directly  observed,  suggests  that  the  upper  limit  is  on  the 
order  of  one  kilometer.  Further  studies  should  be  designed  to  for  these  greater  length  scales. 
The  predetermined  recording  schedule  for  the  CIRCUS  deployment  was  not  serendipitous 
enough  to  acquire  a  teleseismic  event.  Further  studies  should  include  event  triggering  so  that 
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the  gain  for  interesting  seismic  signals  can  be  measured  directly. 

Finally,  the  CIRCUS  array,  being  in  the  horizontal  plane,  made  no  direct  measure¬ 
ments  of  vertical  noise  coherence.  The  present  study  has  successfully  modeled  ambient  noise  as 
horizontally  propagating  trapped  modes,  suggesting  that  these  predominate  over  acoustic 
energy  in  the  water  column.  Future  work  should  test  this  result  more  directly  by  incorporating  a 
vertical  array. 
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Figure  6.1.  CIRCUS  array  resolution  limit  compared  to  expected  seismic  dispersion  curves. 
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Figure  6.2.  Noise  spectrum  while  the  R/\  Melville  and  USNS  Narragansett  are  on  site.  Spec¬ 
tral  lines  are  resolved  with  multiple  taper  harmonic  analysis. 
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Figure  6.3.  Beamfoiming  on  the  ship  signals  at  5.75  Hz.  The  axes  are  wavenumber  times  the 
analysis  frequency  giving  slowness.  Color  at  each  point  is  proportional  to  the  power.  The 
azimuth  of  a  slowness  vector  gives  the  direction  of  propagation  rather  than  the  direction 
towards  the  source.  The  two  squares  are  the  expected  ship  bearings.  At  this  frequency  the  R/V 
Melville  emits  a  strong  signal  and  is  located  successfully. 
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Figure  6.4.  Beamforming  on  the  ship  signals  at  19.5  Hz.  At  this  frequency  the  US  NS  Narra- 
gansett  emits  the  stronger  signal  but  there  is  a  small  error  in  the  bearing  estimation. 
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Figure  6.5.  Beamforming  on  microseism  noise  at  0.25  Hz.  The  observed  phase  slowness  is  not 


resolvably  different  from  zero. 
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Figure  6.6.  Beamforming  on  microseism  noise  at  1.0  Hz.  The  direction  of  propagation  is 
toward  the  southwest  and  with  a  slowness  of  7  s  km-1. 
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Figure  6.7.  Beamforming  on  microseism  noise  at  1.5  Hz.  The  phase  slowness  is  about  1  s  km-1. 
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Figure  6.8.  Summary  of  the  phase  slowness  peaks  for  the  first  20  recording  windows  plotted  as 
a  function  of  frequency.  Scatter  is  substantial  but  the  points  lie  near  the  predicted  dispersion 
curves.  The  solid  part  of  each  curve  is  the  frequency  range  at  which  the  given  mode  is  dom¬ 
inant  according  to  the  coherence  model. 
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Direction  of  Propagation 


A  zimuth 


Figure  6.9.  Summary  of  the  direction  of  microseism  propagation  for  the  same  20  records  as  in 
Figure  6.8.  Azimuth  picks  are  binned  in  10°  intervals.  A  predominant  southwesterly  direction 


is  evident. 
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APPENDIX 

LIST  OF  RECORDED  WINDOWS 


The  stars  in  each  position  indicate  whether  there  is  usable  data  for  the  given  window 
for  each  component  of  each  instrument.  An  empty  circle  indicates  that  there  is  data  for  that 
event  but  that  the  record  start  time  is  significantly  different  from  the  programmed  time.  Note 
that  instruments  1,  2,  4,  8,  12,  13,  and  14  were  programmed  to  record  at  six  hour  intervals 
except  during  the  shooting  phase  which  lasted  from  0200  to  0600  on  day  113,  during  which 
time  they  recorded  every  15  minutes.  Instruments  5  and  6  were  programmed  to  record  at  15 
minute  intervals  until  0600  on  day  1 14  and  every  3  hours  thereafter. 
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**  ** 

2 

112 

00:16 

** 

** 

** 

** 

3 

112 

00:31 

** 

** 

** 

** 

4 

112 

00:46 

** 

** 

** 

** 

5 

112 

01:01 

** 

♦  * 

** 

** 

6 

112 

01:16 

** 

** 

** 

** 

7 

112 

01:31 

** 

*♦ 

** 

** 

8 

112 

01:46 

** 

** 

** 

** 

9 

112 

02:01 

** 

** 

** 

** 

10 

112 

02:16 

** 

** 

** 

** 

11 

112 

02:31 

** 

** 

** 

** 

12 

112 

02:46 

♦* 

** 

** 

** 

13 

112 

03:01 

** 

** 

** 

** 

14 

112 

03:16 

** 

** 

** 

** 

15 

112 

03:31 

** 

** 

** 

** 

16 

112 

03:46 

** 

** 

** 

** 

17 

112 

04:01 

** 

** 

** 

** 

18 

112 

04:16 

** 

** 

** 

** 

19 

112 

04:31 

** 

** 

** 

** 

20 

112 

04:46 

** 

** 

** 

** 

21 

112 

05:01 

** 

** 

** 

** 

22 

112 

05:16  ‘ 

*♦ 

** 

** 

** 

23 

112 

05:31 

** 

** 

** 

** 

24 

112 

05:46 

** 

** 

** 

** 

25 

112 

06:01 

******** 

******** 

******** 

*  ***** 

1 


2 


record 

day 

time  (CUT) 

V 

1  2  4  S  «  S 12  13  14 

P 

1  2  4  J  «  «  12  11 14 

HI 

1  2  4  5  «  S  U  13  14 

H2 

1  2  4  5  «  1 12  13  14 

26 

112 

06:16 

** 

** 

** 

** 

27 

112 

06:31 

** 

** 

** 

** 

28 

112 

06:46 

** 

** 

** 

** 

29 

112 

07:01 

** 

** 

** 

** 

30 

112 

07:16 

** 

** 

** 

** 

31 

112 

07:31 

** 

** 

** 

** 

32 

112 

07:46 

** 

** 

** 

** 

33 

112 

08:01 

** 

** 

** 

** 

34 

112 

08:16 

** 

** 

** 

** 

35 

112 

08:31 

** 

** 

** 

** 

36 

112 

08:46 

** 

** 

** 

** 

37 

112 

09:01 

** 

** 

** 

** 

38 

112 

09:16 

** 

** 

** 

** 

39 

112 

09:31 

** 

** 

** 

** 

40 

112 

09:46 

** 

** 

** 

** 

41 

112 

10:01 

** 

** 

** 

** 

42 

112 

10:16 

** 

** 

** 

** 

43 

112 

10:31 

** 

** 

** 

** 

44 

112 

10:46 

** 

** 

** 

** 

45 

112 

11:01 

** 

** 

** 

** 

46 

112 

11:16 

** 

** 

** 

** 

47 

112 

11:31 

** 

** 

** 

** 

48 

112 

11:46 

** 

** 

** 

** 

49 

112 

12:01 

4c******* 

******* 

******** 

******** 

50 

112 

12:16 

** 

** 

** 

** 

51 

112 

12:31 

** 

** 

** 

** 

52 

112 

12:46 

** 

** 

** 

** 

53 

112 

13:01 

** 

** 

** 

** 

54 

112 

13:16 

** 

** 

** 

** 

55 

112 

13:31 

** 

** 

** 

** 

56 

112 

13:46 

** 

** 

** 

** 

57 

112 

14:01 

** 

** 

** 

** 

58 

112 

14:16 

** 

** 

** 

** 

59 

112 

14:31 

** 

** 

** 

** 

60 

112 

14:46 

** 

** 

** 

** 

61 

112 

15:01 

** 

** 

** 

** 

62 

112 

15:16 

** 

** 

** 

** 

63 

112 

15:31 

** 

** 

** 

** 

64 

112 

15:46 

** 

** 

** 

** 

65 

112 

16:01 

** 

** 

** 

** 

66 

112 

16:16 

** 

** 

** 

** 

67 

112 

16:31 

** 

** 

** 

** 

68 

112 

16:46 

** 

** 

** 

** 

69 

112 

17:01 

** 

** 

** 

** 

70 

112 

17:16 

** 

** 

** 

** 

71 

112 

17:31 

** 

** 

** 

** 

72 

112 

17:46 

** 

** 

** 

** 

73 

112 

18:01 

******** 

******** 

******** 

*  ****** 

74 

112 

18:16 

** 

** 

** 

** 

75 

112 

18:31 

** 

** 

** 

** 

record 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 

121 

122 

123 

124 

125 


3 


day 

time  (CUT) 

V 

1  2  4  5  «  S 12  13  14 

p 

1  2  4  5  «  «  12  1314 

HI 

1  2  4  S  «  S  12  13  14 

H2 

1  2  4  J  «  S  12  13  14 

112 

18:46 

** 

** 

** 

** 

112 

19:01 

** 

** 

** 

** 

112 

19:16 

** 

** 

** 

** 

112 

19:31 

** 

** 

** 

** 

112 

19:46 

%* 

** 

** 

** 

112 

20:01 

** 

** 

** 

** 

112 

20:16 

0* 

0* 

0* 

0* 

112 

20:31 

** 

** 

*★ 

** 

112 

20:46 

** 

** 

** 

** 

112 

21:01 

** 

** 

** 

★* 

112 

21:16 

** 

** 

** 

** 

112 

21:31 

** 

** 

** 

** 

112 

21:46 

** 

** 

** 

** 

112 

22:01 

** 

** 

** 

** 

112 

22:16 

** 

** 

** 

** 

112 

22:31 

** 

** 

** 

** 

112 

22:46 

** 

** 

** 

** 

112 

23:01 

** 

** 

** 

** 

112 

23:16 

★ 

** 

** 

** 

112 

23:31 

** 

** 

** 

** 

112 

23:46 

* 

** 

** 

** 

113 

00:01 

***  **** 

******** 

******** 

*  ****** 

113 

00:16 

** 

** 

** 

113 

00:31 

** 

** 

** 

113 

00:46 

** 

** 

** 

113 

01:01 

** 

** 

** 

113 

01:16 

** 

** 

** 

113 

01:31 

** 

** 

** 

113 

01:46 

** 

** 

** 

113 

02:01 

********* 

********* 

********* 

********* 

113 

02:16 

********* 

********* 

********* 

********* 

113 

02:31 

********* 

********* 

********* 

********* 

113 

02:46 

********* 

********* 

********* 

********* 

113 

03:01 

********* 

********* 

********* 

********* 

113 

03:16 

********* 

********* 

********* 

********* 

113 

03:31 

********* 

********* 

********* 

********* 

113 

03:46 

********* 

********* 

********* 

********* 

113 

04:01 

********* 

********* 

********* 

********* 

113 

04:16 

********* 

********* 

********* 

********* 

113 

04:31 

********* 

********* 

********* 

********* 

113 

04:46 

********* 

********* 

********* 

********* 

113 

05:01 

********* 

********* 

********* 

********* 

113 

05:16 

********* 

********* 

********* 

********* 

113 

05:31 

********* 

********* 

********* 

********* 

113 

05:46 

********* 

********* 

********* 

********* 

113 

06:01 

****  *  ** 

******  ** 

********* 

**  ***  ** 

113 

06:16 

** 

** 

** 

** 

113 

06:31 

** 

** 

** 

** 

113 

06:46 

** 

** 

** 

** 

113 

07:01 

** 

** 

** 

** 

4 


record 

day 

time  (CUT) 

V 

1  2  4  3  <  1 12  1314 

P 

1  243(1121)14 

HI 

1  2  4  J  <  1 12  13  14 

H2 

1  2  4  3  <•  12  13  14 

126 

113 

07:16 

** 

** 

** 

** 

127 

113 

07:31 

** 

** 

** 

** 

128 

113 

07:46 

** 

** 

** 

** 

129 

113 

08:01 

** 

** 

** 

** 

130 

113 

08:16 

** 

** 

** 

** 

131 

113 

08:31 

** 

** 

** 

** 

132 

113 

08:46 

** 

** 

** 

** 

133 

113 

09:01 

** 

** 

** 

** 

134 

113 

09:16 

** 

** 

** 

** 

135 

113 

09:31 

** 

** 

** 

** 

136 

113 

09:46 

** 

** 

** 

** 

137 

113 

10:01 

** 

** 

** 

** 

138 

113 

10:16 

** 

** 

** 

** 

139 

113 

10:31 

** 

** 

** 

** 

140 

113 

10:46 

** 

** 

** 

** 

141 

113 

11:01 

** 

** 

** 

** 

142 

113 

11:16 

** 

** 

** 

** 

143 

113 

11:31 

** 

** 

** 

** 

144 

113 

11:46 

** 

** 

** 

** 

145 

113 

12:01 

**  ****** 

**  ****** 

**  ****** 

******** 

146 

113 

12:16 

** 

** 

** 

** 

147 

113 

12:31 

** 

** 

** 

** 

148 

113 

12:46 

** 

** 

** 

** 

149 

113 

13:01 

** 

** 

** 

** 

150 

113 

13:16 

** 

** 

** 

** 

151 

113 

13:31 

** 

** 

** 

** 

152 

113 

13:46 

★  * 

** 

** 

** 

153 

113 

14:01 

** 

** 

** 

** 

154 

113 

14:16 

** 

** 

** 

** 

155 

113 

14:31 

** 

** 

** 

** 

156 

113 

14:46 

** 

** 

** 

** 

157 

113 

15:01 

** 

** 

** 

** 

158 

113 

15:16 

** 

** 

** 

** 

159 

113 

15:31 

♦ 

** 

** 

** 

160 

113 

15:46 

** 

** 

** 

** 

161 

113 

16:01 

** 

** 

** 

** 

162 

113 

16:16 

** 

** 

** 

** 

163 

113 

16:31 

** 

** 

** 

** 

164 

113 

16:46 

♦  * 

** 

** 

** 

165 

113 

17:01 

** 

** 

** 

** 

166 

113 

17:16 

** 

** 

** 

** 

167 

113 

17:31 

** 

** 

** 

** 

168 

113 

17:46 

** 

** 

** 

** 

169 

113 

18:01 

********* 

********* 

********* 

********* 

170 

113 

18:16 

** 

** 

** 

** 

171 

113 

18:31 

** 

** 

** 

** 

172 

113 

18:46 

** 

** 

** 

** 

173 

113 

19:01 

** 

** 

** 

** 

174 

113 

19:16 

** 

** 

** 

** 

175 

113 

19:31 

** 

** 

** 

** 

record 

176 

177 

178 

179 

180 

181 

182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 

201 

202 

203 

204 

205 

206 

207 

208 

209 

210 

211 

212 

213 

214 

215 

216 

217 

218 

219 

220 

221 

222 

223 

224 

225 


5 


day 

time  (CUT) 

V 

1  1  4  3  1  1121314 

P 

1  2  4  5  S  1 12  1314 

HI 

1  2  4  3  «•  12  13  14 

H2 

1  2  4  3  t  1 12  13  14 

113 

19:46 

** 

** 

** 

** 

113 

20:01 

** 

** 

** 

** 

113 

20:16 

** 

** 

** 

** 

113 

20:31 

** 

** 

** 

** 

113 

20:46 

** 

** 

** 

** 

113 

21:01 

** 

** 

** 

** 

113 

21:16 

** 

** 

** 

** 

113 

21:31 

** 

** 

** 

** 

113 

21:46 

** 

** 

** 

** 

113 

22:01 

** 

** 

** 

** 

113 

22:16 

** 

** 

** 

** 

113 

22:31 

** 

** 

** 

** 

113 

22:46 

** 

** 

** 

** 

113 

23:01 

** 

* 

** 

** 

113 

23:16 

** 

** 

** 

** 

113 

23:31 

** 

** 

** 

** 

113 

23:46 

** 

** 

** 

** 

114 

00:01 

********* 

********* 

********* 

********* 

114 

00:16 

** 

** 

** 

** 

114 

00:31 

** 

** 

** 

** 

114 

00:46 

** 

** 

** 

** 

114 

01:01 

*★ 

** 

** 

** 

114 

01:16 

** 

** 

** 

** 

114 

01:31 

** 

** 

** 

** 

114 

01:46 

** 

** 

** 

** 

114 

02:01 

** 

** 

** 

** 

114 

02:16 

** 

** 

** 

** 

114 

02:31 

** 

** 

** 

** 

114 

02:46 

** 

* 

** 

** 

114 

03:01 

** 

* 

** 

** 

114 

03:16 

** 

* 

** 

** 

114 

03:31 

** 

* 

** 

** 

114 

03:46 

** 

* 

** 

** 

114 

04:01 

** 

* 

** 

** 

114 

04:16 

** 

* 

** 

** 

114 

04:31 

** 

* 

** 

** 

114 

04:46 

** 

* 

** 

** 

114 

05:01 

** 

* 

** 

** 

114 

05:16 

** 

* 

** 

** 

114 

05:31 

** 

* 

** 

** 

114 

05:46 

* 

* 

** 

** 

114 

06:01 

******  ** 

***  ***** 

******  ** 

******  ** 

114 

09:01 

** 

* 

** 

** 

114 

12:01 

**  ***** 

**  ***** 

** 

****** 

** 

****** 

114 

15:01 

* 

* 

** 

** 

114 

18:01 

**  **Q*Q* 

**  *0*0* 

** 

**o*o* 

** 

**0*0* 

114 

21:01 

** 

* 

** 

** 

115 

00:01 

**  ***  ** 

**  **  ** 

** 

****** 

** 

***  ** 

115 

03:01 

** 

* 

** 

* 

115 

06:01 

**  ***** 

**  ***** 

** 

0***** 

** 

0***** 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 

261 

262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 

273 

274 

275 


day 

time  (CUT) 

V 

i  a  4  s  •  •  i*  n  M 

p 

1  1  4  5  «  •  11 13  14 

HI 

1  2  4  S  «  1 12  1314 

H2 

I  2  4  5  4  B  12  13  14 

115 

09:01 

** 

* 

** 

** 

115 

12:01 

**  ****** 

** 

***** 

**  ****** 

**  ****** 

115 

15:01 

** 

* 

** 

** 

115 

18:01 

**  ***** 

** 

***** 

**  ****** 

**  ****** 

115 

21:01 

+* 

* 

** 

** 

116 

00:01 

*Q  ***  ** 

*0 

**  ** 

*Q  ***  ** 

*0  ***  ** 

116 

03:01 

** 

* 

** 

** 

116 

06:01 

*****Q*** 

** 

0*** 

****  Q*** 

****  q*** 

116 

09:01 

** 

* 

** 

** 

116 

12:01 

**Q****** 

**0 

***** 

**Q****** 

**0****** 

116 

15:01 

** 

* 

** 

** 

116 

18:01 

***  ***** 

*** 

***** 

********* 

********* 

116 

21:01 

** 

* 

** 

** 

117 

00:01 

*******  * 

*** 

***** 

********* 

********* 

117 

03:01 

** 

★ 

** 

** 

117 

06:01 

**Q**  *+* 

**0 

*  *** 

**Q**  *** 

**0**  *** 

117 

09:01 

** 

★ 

** 

** 

117 

12:01 

********* 

*** 

***** 

********* 

********* 

117 

15:01 

** 

* 

** 

** 

117 

18:01 

********* 

*** 

***** 

***  ***** 

***  ***** 

117 

21:01 

* 

* 

0* 

0* 

118 

00:01 

***  **  ** 

*** 

**  ** 

******  ** 

***  **  ** 

118 

03:01 

** 

* 

** 

* 

118 

06:01 

**Q****** 

**0 

**  ** 

**Q***  ** 

**0***  ** 

118 

09:01 

** 

* 

** 

** 

118 

12:01 

** 

**  ** 

*Q****  ** 

*0****  ** 

118 

15:01 

** 

* 

** 

** 

118 

18:01 

**Q***  ** 

**0 

*+  %* 

**Q***  ** 

**0***  ** 

118 

21:01 

** 

* 

** 

** 

119 

00:01 

********Q 

*** 

**** 

********0 

********0 

119 

03:01 

** 

* 

** 

** 

119 

06:01 

Q******** 

0** 

***  * 

Q******  * 

0******  * 

119 

09:01 

** 

* 

** 

** 

119 

12:01 

******  ** 

*** 

**  ** 

********* 

******  ** 

119 

15:01 

** 

♦ 

** 

** 

119 

18:01 

Q*****  ** 

**  ** 

Q******** 

0*****  ** 

119 

21:01 

** 

♦ 

** 

** 

120 

00:01 

********* 

*** 

***** 

********* 

********* 

120 

03:01 

** 

% 

** 

** 

120 

06:01 

******  ** 

*** 

**  ** 

***  **  ** 

***  **  ** 

120 

09:01 

** 

% 

** 

** 

120 

12:01 

********* 

*** 

***** 

********* 

********* 

120 

15:01 

** 

* 

** 

** 

120 

18:01 

********Q 

*** 

****Q 

********Q 

********0 

120 

21:01 

0* 

* 

0 * 

** 

121 

00:01 

*****Q*** 

*** 

+0++* 

*****0*** 

*****0*** 

121 

03:01 

** 

* 

** 

** 

121 

06:01 

********* 

*** 

***** 

********* 

********* 

121 

09:01 

*0 

0 

0 

*0 

121 

12:01 

********* 

*** 

***** 

********* 

********* 

record 

276 

277 

278 

279 

280 

281 

282 

283 

284 

285 

286 

287 

288 

289 

290 

291 

292 

293 

294 

295 

296 

297 

298 

299 

300 

301 

302 

303 

304 

305 

306 

307 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 

321 

322 

323 

324 

325 


7 


day 

time  (CUT) 

V 

1  2  4  J  4  4  12  13  14 

P 

1  2  4  5  C  1 12  13 14 

HI 

1  2  4  5  4  4 12  13  14 

H2 

1  2  4  5  4  4  12  13  14 

121 

15:01 

** 

* 

** 

** 

121 

18:01 

******  ** 

*** 

***** 

********* 

********* 

121 

21:01 

** 

* 

** 

** 

122 

00:01 

********* 

*** 

***** 

********* 

********* 

122 

03:01 

** 

* 

** 

** 

122 

06:01 

********* 

*** 

**  ** 

********* 

********* 

122 

09:01 

** 

* 

* 

** 

122 

12:01 

****Q**** 

*** 

Q**** 

****Q**** 

****Q**** 

122 

15:01 

** 

* 

★ 

** 

122 

18:01 

*0**0**** 

*0* 

Q**** 

*Q**  **** 

*0**0**** 

122 

21:01 

** 

* 

* 

** 

123 

00:01 

******  ** 

*** 

***** 

********* 

******  ** 

123 

03:01 

** 

* 

** 

** 

123 

06:01 

******  ** 

*** 

**  ** 

******  ** 

******  ** 

123 

09:01 

** 

* 

** 

** 

123 

12:01 

Q******** 

0** 

***** 

Q******** 

Q******** 

123 

15:01 

** 

* 

ft* 

** 

123 

123 

18:01 

21:01 

*****Q**Q 

0 

*** 

*0  *0 

*****Q**Q 

0 

*****Q**0 

*0 

124 

00:01 

********Q 

*** 

****Q 

********Q 

********Q 

124 

03:01 

** 

* 

** 

** 

124 

06:01 

*Q****  ** 

*0* 

**  ** 

*Q****  ** 

*Q****  ** 

124 

09:01 

ft* 

* 

** 

ft* 

124 

12:01 

********* 

*** 

***** 

********* 

********* 

124 

15:01 

** 

* 

** 

** 

124 

18:01 

******  ** 

*** 

***** 

********* 

********* 

124 

21:01 

** 

* 

♦  ft 

ft* 

125 

00:01 

********* 

*** 

***** 

********* 

********* 

125 

03:01 

** 

* 

** 

** 

125 

06:01 

********* 

*** 

***** 

********* 

********* 

125 

09:01 

ft* 

* 

** 

** 

125 

12:01 

*******Q* 

*** 

***Q* 

*******Q* 

*******Q* 

125 

15:01 

** 

* 

It* 

ft* 

125 

18:01 

********* 

♦** 

***** 

********* 

********* 

125 

21:01 

** 

* 

** 

** 

126 

00:01 

*Q****  ** 

*0* 

***** 

*Q****  ** 

*Q******* 

126 

03:01 

** 

* 

ft* 

** 

126 

06:01 

*****Q  ** 

*** 

*0  ** 

*****Q  ** 

*****Q  ** 

126 

09:01 

** 

* 

* 

** 

126 

12:01 

*****Q*** 

*** 

*Q*** 

*****Q*** 

*****Q*** 

126 

15:01 

* 

* 

** 

** 

126 

18:01 

********* 

*** 

***** 

********* 

********* 

126 

21:01 

♦ 

* 

ft* 

★  * 

127 

00:01 

****  **** 

★** 

***** 

********* 

********* 

127 

03:01 

* 

* 

** 

** 

127 

06:01 

***  *  ** 

*** 

**  ** 

***Q**  ** 

***Q**  ** 

127 

09:01 

* 

* 

** 

** 

127 

12:01 

***  *  ** 

*** 

**  ** 

******  ** 

******  ** 

127 

15:01 

★ 

* 

0* 

0* 

127 

18:01 

****  **** 

*** 

***  * 

********* 

*******  * 
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record 

day 

time  (CUT) 

V 

1  2  4  5  4  •  12  13  14 

p 

1  2  4  J  «  1 12  13  M 

HI 

1  2  4  3  *  1 12  13  14 

H2 

1  2  4  3  «  1 12  13  14 

326 

127 

21:01 

* 

* 

** 

** 

327 

128 

00:01 

**** 

* 

** 

*** 

**  ** 

******  ** 

******  ** 

328 

128 

03:01 

* 

* 

* 

329 

128 

06:01 

**0* 

* 

** 

** 

*  ** 

**0* 

**** 

**o* 

*  ** 

330 

128 

09:01 

* 

* 

* 

331 

128 

12:01 

**** 

* 

** 

*** 

**** 

**** 

**** 

***% 

***% 

332 

128 

15:01 

* 

* 

* 

333 

128 

18:01 

*%%* 

**** 

*** 

**** 

**  * 

**** 

*%* 

*  ** 

334 

128 

21:01 

* 

* 

* 

335 

129 

00:01 

*%** 

* 

*0 

*** 

*  *0 

**** 

*  *0 

*%** 

*  *0 

336 

129 

03:01 

* 

* 

* 

337 

129 

06:01 

**** 

* 

** 

*** 

*  ** 

**** 

*  ** 

***% 

*  ** 

338 

129 

09:01 

* 

* 

* 

339 

129 

12:01 

**** 

* 

** 

** 

*  ** 

**  * 

*  ** 

**  * 

*  *% 

340 

129 

15:01 

* 

* 

* 

341 

129 

18:01 

**** 

**** 

*** 

**  * 

**** 

**** 

*** 

**** 

342 

129 

21:01 

* 

343 

130 

00:01 

**** 

* 

00 

*** 

**00 

**** 

*  0 

**** 

*  0 

344 

130 

03:01 

* 

* 

* 

345 

130 

06:01 

**** 

***o 

*** 

*0 

**%* 

*0 

**% 

*0 

346 

130 

09:01 

* 

* 

* 

347 

130 

12:01 

**** 

**** 

*** 

**** 

**** 

**** 

**** 

348 

130 

15:01 

* 

* 

* 

349 

130 

18:01 

**** 

**** 

*** 

**** 

**** 

**** 

**** 

**** 

350 

130 

21:01 

♦ 

* 

351 

131 

00:01 

**** 

* 

** 

*** 

*  *% 

*** 

**** 

*** 

★  ** 

352 

131 

03:01 

* 

353 

131 

06:01 

**** 

* 

** 

*** 

*  ** 

*** 

*** 

*** 

*  ** 

354 

131 

09:01 

* 

* 

355 

131 

12:01 

*0*0 

* 

0* 

*0* 

*  * 

*0* 

*  0 

*0*0 

*  0* 

356 

131 

15:01 

* 

* 

357 

131 

18:01 

%%** 

0 

%* 

%*% 

Q#  +  * 

** 

0*** 

*** 

0  * 

358 

131 

21:01 

* 

359 

132 

00:01 

**** 

0 

** 

%*% 

0  ** 

%%* 

0  ** 

*** 

0  ** 

360 

132 

03:01 

* 

361 

132 

06:01 

*0  * 

* 

** 

*0 

*  ** 

*0 

*  * 

*0 

*  ** 

362 

132 

09:01 

* 

363 

132 

12:01 

**** 

* 

*% 

*** 

*  ** 

**% 

*  ** 

*** 

*  ** 

364 

132 

15:01 

★ 

365 

132 

18:01 

***% 

* 

%% 

*%* 

** 

%** 

** 

*** 

** 

366 

132 

21:01 

* 

* 

367 

133 

00:01 

%*%* 

* 

** 

*** 

*  * 

*** 

*  ** 

**** 

*  ** 

368 

133 

03:01 

* 

369 

133 

06:01 

**** 

* 

** 

%** 

*  * 

%%* 

*  ** 

*** 

*  ** 

370 

133 

09:01 

* 

371 

133 

12:01 

**** 

* 

0* 

*** 

*  * 

%** 

*  o* 

**** 

*  o* 

372 

133 

15:01 

* 

* 

373 

133 

18:01 

**** 

* 

** 

*** 

*  * 

*** 

*  ** 

**** 

*  ** 

374 

133 

21:01 

* 

375 

134 

00:01 

**** 

* 

** 

*** 

*  * 

*** 

*  ** 

*** 

*  ** 

record 

376 

377 

378 

379 

380 

381 

382 

383 

384 

385 

386 

387 

388 

389 

390 

391 

392 

393 

394 

395 

396 

397 

398 

399 

400 

401 

402 

403 

404 

405 

406 

407 

408 

409 

410 

411 

412 

413 

414 

415 

416 

417 

418 

419 

420 

421 

422 

423 

424 

425 
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day 

time  (CUT) 

V 

124  JIIU»14 

P 

1  2  4  J  ••  12  13  14 

HI 

1  2  4  J  t  1 12  13  14 

H2 

1  2  4  3  S  S  12  13  14 

134 

134 

03:01 

06:01 

* 

**** 

**** 

4c  4c 

4c  4c  4c 

4c  4c 

4*  4*  4*  4* 

4*4* 

4*  4*  4*  4* 

134 

134 

09:01 

12:01 

0 

**** 

*  ** 

4c  4c 

4c  4c 

4c  4c  4c 

4*  4*4* 

4*4*4* 

4*  4*4* 

134 

134 

15:01 

18:01 

* 

*0** 

4*  4*4* 

*0* 

4c  4c 

*0* 

4*  4*4* 

4*0* 

4*  4* 

134 

135 

21:01 

00:01 

* 

4*  4*  4*  4* 

*  ** 

4c  4c  4c 

4c  4c 

4c  4c  4c 

4*  4* 

4*4*4* 

4*  4*4* 

135 

135 

03:01 

06:01 

* 

*** 

*  ** 

4c  4c  4c 

4c  4c 

4c  4c  4c 

4*  4*4* 

4*4*4* 

4*  4*4* 

135 

135 

09:01 

12:01 

*** 

*  ** 

4c  4c  4c 

4c  4c 

4c  4c  4c 

4*  4*4* 

4*4* 

4*  4*4* 

135 

135 

15:01 

18:01 

*0* 

♦  0* 

4cq 

4c 

4*o 

*  0* 

4*0 

*  0* 

135 

136 

21:01 

00:01 

0** 

*  4*4* 

o* 

4c  4c 

0* 

4*  4*4* 

o** 

4*  4*4* 

136 

136 

03:01 

06:01 

*0* 

4*4*4*4* 

*0 

4c  4c  4c 

*0 

4*  4*  4*  4* 

4*o 

4*  4*  4*  4* 

136 

136 

09:01 

12:01 

** 

*  ** 

4c  4c 

4c  4c 

4*4* 

4*  4*4* 

4*4* 

4*  4*4* 

136 

136 

15:01 

18:01 

** 

**** 

4c  4c 

4c  4c  4c 

4*4* 

4*  4*  4*  4* 

4*4* 

4*  4*  4*  4* 

136 

137 

21:01 

00:01 

** 

*  ** 

4c  4c 

4c  4c 

4*4* 

4*  4*  4*  4* 

4*4* 

4*  4*4* 

137 

137 

03:01 

06:01 

** 

*  ** 

4c  4c 

4c  4c 

4*4* 

4*  4*4* 

4*4* 

4*  4*4* 

137 

137 

09:01 

12:01 

** 

*  ** 

4c  4c 

4c  4c 

4*4* 

4* 

4* 

4*  4* 

137 

137 

15:01 

18:01 

** 

*  ** 

4c  4c 

4c 

4*4* 

4*  4* 

4*4* 

4*  4* 

137 

138 

21:01 

00:01 

** 

*  ** 

4c  4c 

4c  4c 

4*4* 

4*  4*4* 

4*4* 

4*  4*4* 

138 

138 

03:01 

06:01 

** 

4*  ** 

4c  4c 

4c  4c 

4*4* 

4*  4*4* 

4*4* 

4*  4*4* 

138 

138 

09:01 

12:01 

** 

♦  ** 

4c  4c 

4c  4c 

4*4* 

4*4*  4* 

4*4* 

4*  4*4* 

138 

138 

15:01 

18:01 

** 

4*  4c 

4c  4c 

4*4* 

4* 

4*4* 

4* 

138 

139 

21:01 

00:01 

** 

4c  4c 

4c  4c 

4c 

4*4* 

4*  4* 

4*4* 

4*  4* 

139 

139 

03:01 

06:01 

4*  4* 

4c4c  4c 

4c  4c 

4c  4c 

4*4* 

4*4*  4* 

4*4* 

4*4*  4* 

139 

139 

09:01 

12:01 

** 

*  0 

4c  4c 

4c 

4*4* 

4*  o 

4*4* 

4*  o 

139 

139 

15:01 

18:01 

4*  4* 

4c  4c 

4c  4c 

4c 

4*4* 

4*  4* 

4*4* 

4*  4* 

139 

140 

21:01 

00:01 

4*  4* 

4*4*  4c 

4c  4c 

4c 

4*4* 

4*  4* 

4*4* 

4*  4* 

140 

140 

03:01 

06:01 

** 

4c  4c 

4c  4c 

4c 

4*4* 

4*  4* 

4*4* 

4*  4* 

10 


record 


426 

427 

428 

429 

430 

431 

432 

433 

434 

435 

436 

437 

438 

439 

440 

441 

442 

443 

444 

445 

446 

447 

448 

449 

450 

451 


day 

time  (CUT) 

I  2  4 

V 

3  *  «  12 13  14 

P 

1  2  4  S  *  «  12  13  14 

HI 

1  2*  311111314 

1  2  4 

H2 

5  •  1 12  1314 

140 

140 

09:01 

12:01 

*0 

4*  4* 

*0 

4* 

*0 

4*  4* 

*0 

4*4* 

4* 

140 

140 

15:01 

18:01 

4*4* 

4*  4* 

4*4* 

4* 

4*4* 

4* 

4*4* 

4* 

4* 

140 

141 

21:01 

00:01 

*0 

4*  4* 

*0 

4* 

*0 

4*4*  4* 

*0 

4* 

4* 

141 

141 

03:01 

06:01 

0* 

4*  4* 

0* 

4*4* 

0* 

4*4*  4* 

0* 

4*4* 

4* 

141 

141 

09:01 

12:01 

4*4* 

4*4*  4* 

4*4* 

4*4* 

4*4* 

4*4* 

4*4* 

4*4* 

4* 

141 

141 

15:01 

18:01 

4*4* 

4* 

4*4* 

4* 

4*4* 

4* 

4*4* 

4* 

141 

142 

21:01 

00:01 

4*4* 

0 

4*4* 

0* 

4*4* 

0* 

4*4* 

0* 

142 

142 

03:01 

06:01 

4*4* 

4*4* 

4*4* 

4*4* 

4*4* 

4*4* 

4*4* 

4*4* 

142 

142 

09:01 

12:01 

0* 

4*4* 

0* 

4*4* 

0* 

4*4* 

0* 

4*4* 

142 

142 

15:01 

18:01 

4*4* 

4* 

4*4* 

4* 

4*4* 

4* 

4*4* 

4* 

142 

143 

21:01 

00:01 

4*4* 

4*4* 

4*4* 

4*4* 

4*4* 

4*4* 

4*4* 

4*4* 

143 

143 

03:01 

06:01 

4*4* 

4*4* 

4*4* 

4*4* 

4* 

4* 

4*4* 

4*4* 

143 

143 

09:01 

12:01 

4*4* 

4*4* 

4* 

4*4* 
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